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ABSTRACT 

Recent  experimental  studies  strongly  ^suggest  that  the  hairpin 
vortex  plays  an  important,  a«d  possibly  dominant,  role  in  the  dynamics 
of  turbulent  flows  near  walls.  In  this  study,  various  aspects  of  the 
motion  of  hairpin  vortices  near  solid  walls  are  addressed  on  a  theo¬ 
retical  basis.  In  particular  the  following  topics  are  of  interest: 
(1)  the  nature  of  the  evolution  of  hairpin  vortices  in  a  shear  flow  ; 
{2y  the  type  of  flow  induced  near  a  wall  by  a  convected  hairpin 
vortex;  (3)  the  character  of  the  viscous  response  near  a  wall  to  the 
hairpin  vortex  motion  and  (4)  the  nature  of  the  interaction  of  two  ' 
hairpin  vortices. 

In  the  first  phase  of  this  study,  a  numerical  procedure  is 
developed  to  allow  the  accurate  evaluation  of  the  trajectory  of  a 
three-dimensional^  vortex  for  vortices  having  small  cores.  The 
integration  method  is  based  on  a  numerical  approximation  to  the 
Biot-Savart  integral;  most  existing  vortex  calculation  methods  have 
severe  stability  problems  for  vortices  with  small  cores.  The- 
stability  problem  fs  overcome  with  the  present  method  and  the  tech¬ 
nique  is  applied  to  compute  the  evolution  of  convected  vortex  loops 
and  hairpin  vortices,  both  In  a  uniform  flow  and  in  a  shear  flow 
above  a  wall.  For  the  case  of*> hairpin  vortices  evolving  in  a  shear 
flow,  a  regenerative  process  Is  observed  wherein  secondary  hairpin 
vortices  form  outboard  of  the  original  hairpin  vortex  in  a  manner 
consistent  with  experimental  observations.  Calculations  for  two 
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hairpin  vortices,  which  are  Initially  close  to  one  another  in  a 
shear  flow,  show  that  the  vortices  interact  In  a  manner  so  as  to 
reinforce  one  another.  Simulations  are  carried  out  to  compute  the 
evolution  of  hydrogen  bubble  lines  as  a  hairpin  vortex  passes  and 
these  show  that  such  vortices  can  give  rise  to  wall  layer  streaks. 

In  the  second  phase  of  this  study,  the  response  is  computed 
for  a  viscous  flow  near  a  wall  due  to  the  motion  of  a  hairpin 
vortex  above  the  wall.  The  results  reveal  that  a  complex,  unsteady 
boundary* layer  flow  develops  near  the  wall  which  ultimately  develops 
a  very  strong  local  outward  growth.  The  final  stages  of  this 
development  are  expected  to  lead  to  a  boundary- layer  eruption  and 
the  creation  of  a  secondary  hairpin  vortex  by  a  strong  viscous- 
invlscid  interaction  with  the  outer  flow.  This  mode  of  regeneration 
is  believed  to  be  a  fundamental  process  wherein  new  vorticity  from 
the  wall  region  is  continually  introduced  into  the  turbulent  bound¬ 
ary  layer. 
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Figure  3.12 

-  End  view  of  the  development  of 
simulated  hydrogen  bubble  lines  due 
to  a  convected  ring  (V=3). 
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Figure  3.13 

-  Side  view  of  the  trajectory  of  a 
vortex  loop  convected  with  speed 

V=10  (a=135°). 
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Figure  3.14 

-  Temporal  development  of  the  simulated 
hydrogen  bubble  lines  due  to  a 
convected  vortex  loop  (V=10). 
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(c)  t=l .  and  (d)  t=2. 
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Figure  3.15 

-  Side  view  of  the  development  of 
simulated  hydrogen  bubble  lines  due 
to  a  convected  vortex  loop  with  speed 

V=10  (a=l 35°). 
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Figure  3.16 

-  End  view  of  the  development  of 
simulated  hydrogen  bubble  lines  due 
to  a  convected  vortex  loop  with  speed 

V=10  (a=l  35°). 
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Figure  4.1 

-  Vortex  motion  above  a  wall. 
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-  Temporal  development  for  a  hairpin 
vortex  In  a  stagnant  flow.  Case  1 
(a=45°,  A=0.5,  8=20);  the  vortex 
position  is  plotted  every  40  time 
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Figure  4.3 

-  Temporal  development  for  a  hairpin 
vortex  in  a  stagnant  flow.  Case  2 
(a=45° ,  A=l.,  8=20);  the  vortex 
position  is  plotted  every  20  time 
steps  ( At=0 . 0002). 
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Figure  4.4 

-  Temporal  development  for  a  hairpin 
vortex  In  a  stagnant  flow.  Case  3 
(a*0°,  A=0.5,  8=20);  the  vortex 
position  is  plotted  every  30  time 
steps  (At=0.0002). 
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Temporal  development  for  a  hairpin 
vortex  in  a  shear  flow  (V=50).  Case  4 
(a=45° ,  A=0.5,  6=20);  the  vortex 
position  is  plotted  every  20  time 
steps  (At=0.0005). 

Temporal  development  for  a  hairpin 
vortex  in  a  shear  flow  (V=150).  Case  5 
(a=45°,  A=0.5,  6=20);  the  vortex 
position  Is  plotted  every  20  time 
steps  (At=0.0003). 

Temporal  development  for  a  hairpin 
vortex  in  a  shear  flow  (V=250).  Case  6 
(a=45° ,  A=0.5,  6=20);  the  vortex 
position  is  plotted  every  30  time 
steps  (At=0.0002). 

Continuation  of  Integration  for  Case  6; 
the  vortex  position  is  plotted  every 
15  time  steps  (At=0.0002). 

Flow  visual Izatlon  sequence  reproduced 
from  Acalar  and  Smith  (1984)  showing  a 
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The  simulated  hydrogen  bubble  lines  due 
to  a  convected  hairpin  vortex  (Case  6) 
at  t  =0.3;  initial  bubble  lines  were 
located  at  y  =  .8. 

The  simulated  hydrogen  bubble  lines  due 
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at  t  =  .03;  initial  bubble  lines  were 
located  at  y  =  .7. 
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Figure  4.19(a)  - 


The  simulated  hydrogen  bubble  lines 
due  to  a  convected  hairpin  vortex 
(Case  6)  at  t  =  .06;  initial  bubble 
lines  were  located  at  y  =  .7. 

The  simulated  hydrogen  bubble  lines 
due  to  a  convected  hairpin  vortex 
(Case  6)  at  t  =  .07;  initial  bubble 
lines  were  located  at  y  =  .7. 

Side  view  of  hairpin  vortex  (at 
t  =  0.1  for  Case  4  of  Section  4.5) 
and  instantaneous  streamwise  measur¬ 
ing  locations  . 

Instantaneous  streamwise  velocity 
fluctuations  at  stations  Indicated 
in  Figure  4.16(a). 

Instantaneous  streamwise  velocity 
profiles  at  stations  Indicated 
in  Figure  4.16(a). 

Side  view  of  the  developing  hairpin 
vortex  for  Case  4  at  various  times. 

Instantaneous  velocity  fluctuations 
of  the  measuring  station  indicated 
in  Figure  4.17(a);  note  the  shifted 
origins. 

Instantaneous  velocity  profiles  at 
the  measuring  station  indicated  in 
Figure  4.17(a). 

Side  view  of  hairpin  vortex  (at 
t  »  0.78  for  Case  6  of  Section  4.5) 
and  instantaneous  streamwise  measur¬ 
ing  locations. 

Instantaneous  streamwise  velocity 
fluctuations  at  stations  indicated 
In  Figure  4.18(a). 

Side  view  of  the  developing  hairpin 
vortex  for  Case  6  at  various  time- 
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4.19(b) 

-  Instantaneous  velocity  fluctuations 
at  the  measuring  station  indicated 
in  Figure  4.19(a);  note  the  shifted 
origins- 
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fluctuation  of  hairpin  vortex  over  a 
cycle  as  obtained  from  hydrogen  bubble 
flow  visualization  pictures  on  plane  of 
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-  Development  of  the  two  hairpin 
vortices  at  t  =  .015. 
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-  Development  of  the  two  hairpin 
vortices  at  t  =  .04. 
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-  Development  of  the  two  hairpin 
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Figure  5.6 


Temporal  evolution  of  e,  and  e2>  the 
decomposition  of  0  . 


Figure  5.7(a)  - 


Streamlines  In  the  symmetry  plane  at 
t  =  0.001. 


Figure  5.7(b)  -  Lines  of  constant  e  at  t  =  0.001. 


Figure  5.7(c) 


Limiting  surface  streamlines  in  the 
xz  plane  at  t  =  0.001. 


Figure  5.8(a)  - 


Streamlines  In  the  syrrmetry  plane  at 
t  =  0.01. 


Figure  5.8(b)  -  Lines  of  constant  e  at  t  =  0.01. 


Figure  5.8(c) 


Limiting  surface  streamlines  In  the 
xz  plane  at  t  =  0.01. 


Figure  5.9(a) 


Streamlines  In  the  symmetry  plant  at 
t  =  0.02. 


Figure  5.9(b)  - 

Figure  5.9(c)  - 


Lines  of  constant  e  at  t  =  0.02. 


Limiting  surface  streamlines  In  the 
xz  plane  at  t  -  Q.02. 


Figure  5.10(a)  - 


Streamlines  in  the  syirmetry  plane  at 
t  =  0.03. 


Figure  5.10(b)  -  Lines  of  constant  0  at  t  =  0.03 


Figure  5.10(c)  - 


Limiting  surface  streamlines  in  the 
xz  plane  at  t  *  0. 03- 


Figure  5.11(a) 


Streamlines  in  the  symmetry  plane  at 
t  =  0.04. 


Figure  5.11(b)  -  Lines  of  constant  e  at  t  =  0.04. 


Figure  5.11(c)  - 


Limiting  surface  streamlines  in  the 
xz  plane  at  t  =  0.04. 


Figure  5.12(a)  - 


Streamlines  in  the  symmetry  plane  at 
t  =  0.05. 


Figure  5.12(b)  -  Lines  of  constant  0  at  t  =  0.05. 


Figure  5.12(c) 


Limiting  surface  streamlines  in  the 
xz  plane  at  t  =  0.05. 
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Figure  5.13(a)  -  Streamlines  in  the  symmetry  plane  at 
t  =  0.06. 

Figure  5.13(b)  -  Lines  of  constant  e  at  t  =  0.06. 

Figure  5.13(c)  -  Limiting  surface  streamlines  in  the 
xz  plane  at  t  =  0.06. 

Figure  5.14(a)  -  Streamlines  in  the  symmetry  plane  at 
t  =  0.07. 

Figure  5.14(b)  -  Lines  of  constant  e  at  t  =  0.07. 

Figure  5.14(c)  -  Limiting  surface  streamlines  in  the 
xz  plane  at  t  =  0.07. 


Figure  5.15 


Temporal  development  of  displacement 
thickness. 


CHAPTER  I 


INTRODUCTION 

1.1  Background 

Most  flows  occurring  in  nature  and  in  engineering  application 
are  turbulent.  The  atmospheric  boundary  layer  as  well  as  ocean 
currents  are  generally  in  turbulent  motion.  The  major  portion  of 
boundary  layers  which  form  and  grow  on  such  devices  as  turbine  blades 
and  aircraft  wings  are  turbulent.  Turbulence  is  usually  present  in 
most  channel  and  pipe  flows  as  well  as  In  most  combustion  processes. 
The  study  and  prediction  of  turbulent  flow  phenomena  is  thus  an 
extremely  important  area  of  basic  research. 

Generally  turbulence  cannot  sustain  itself  without  some  source 
of  energy.  For  incompressible  flows,  the  production  of  new  turbu¬ 
lence  Is  generally  observed  to  be  associated  with  the  presence  of 
shear  in  the  mean  flow;  the  most  common  location  where  mean  shear 
occurs  is  near  solid  walls  where  the  flow  must  be  reduced  to  rela¬ 
tive  rest  to  satisfy  the  no  slip  condition.  The  role  of  shear  in 
the  production  process  appears  to  be  important  for  if  turbulence 
arrives  in  an  environment  where  there  is  no  shear  (or  some  other 
maintenance  mechanism),  it  decays.  Unfortunately,  the  production 
process  is  not  well  understood.  A  common  view  is  that  vortex 
motions  play  an  important  role  in  the  transport  of  turbulence;  the 
large  scale  eddies  are  thought  to  extract  energy  from  the  mean  flow 


(see  for  example  Tennekes  and  Lumley,  1972,  p.  41  or  Hlnze,  1975, 
p.  74)  and  then  to  transfer  energy  to  smaller  scale  eddies  through 
vortex  stretching.  Ultimately  the  smaller  eddies  dissipate  energy 
Into  heat  through  viscous  dissipation.  Although  this  statement  of 
the  processes  In  turbulence  Is  commonly  accepted.  It  does  not  shed 
much  light  on  the  complex  dynamics  that  occur  in  the  production 
mechanisms  near  walls. 

In  recent  years,  there  has  been  considerable  experimental 
Interest  In  the  structure  of  turbulent  boundary -layer  flow  near 
walls  and  it  has  become  apparent  that  such  flows  are  not  comprised 
of  random  fluctuations  in  velocity;  rather  an  ordered  (but  complex) 
cyclic  sequence  of  events  appears  to  take  place  in  the  flow  near 
the  wall.  This  apparently  deterministic  sequence  Is  usually 
referred  to  as  the  coherent  structure  or  the  coherent  behavior  of 
the  wall-layer  flow;  here  the  wall  layer  Is  understood  to  imply  that 
portion  of  a  turbulent  flow  between  the  wall  and  locations  where  the 
mean  profile  Is  logarithmic  In  y,  the  distance  from  the  wall.  There 
are  two  features  which  are  observed  to  dominate  the  dynamics  of  the 
wall- layer  flow,  namely  the  low  speed  streaks  and  the  bursting 
phenomenon.  When  a  visualization  medium,  such  as  dye  or  hydrogen 
bubbles,  is  introduced  into  the  wall -layer  flow  it  is  observed  to 
collect  into  low-speed  streaks  near  the  wall  which  are  aligned  in 
the  streamwise  direction  and  which  persist  over  relatively  long 
distances;  the  streaks  are  referred  to  as  "low-speed"  because  the 
streamwise  velocity  near  the  streaks  is  generally  observed  to  be 


less  than  the  velocity  distribution  corresponding  to  the  mean  pro¬ 
file.  A  schematic  sketch  of  the  observed  phenomenon  is  given  in 
Figure  1.1. 

The  second  important  feature  of  the  wall -layer  flow  is  the 
bursting  process  wherein  the  wall -layer  fluid  is  observed  to  erupt, 
at  Isolated  spanwlse  and  streanwlse  locations,  into  the  outer  region 
of  the  flow.  The  event  is  observed  to  initiate  near  a  wall  layer 
streak  and  in  its  primary  stages  is  characterized  by  a  lifting  of 
the  low-speed  streak  from  the  wall;  an  apparently  oscillatory 
behavior  then  starts  to  develop  In  the  streak  and  then  a  violent 
chaotic  eruption  of  wall-layer  fluid  is  observed.  A  schematic 
diagram  of  the  various  stages  of  the  process,  as  viewed  from  the 
side  and  which  is  adapted  from  Kline  (1978),  is  shown  in  Figure 
1.2. 

The  physical  cause  of  the  wall-layer  streaks  and  the  bursting 
phenomenon  is  not  well  understood.  In  this  study,  the  possibility 
that  the  observed  effects  are  due  to  convected  vortex  motion  will  be 
considered.  In  particular  the  major  thrust  of  this  study  is  an 
investigation  of  the  cause  and  effect  relationships  generated  by  a 
convecting  hairpin  vortex. 

1.2  Experimental  Studies 

It  is  commonly  believed  that  the  vortex  motions  that  are 
clearly  observable  in  boundary-layer  turbulence,  play  an  important 
role  in  the  dynamics  of  the  turbulence  and  a  variety  of  model  vortex 


Figure  1.1  -  A  schematic  sketch  of  dye  collected 
into  low-speed  streaks  in  a  turbulent 
flow  near  a  wall. 
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motions  have  been  suggested.  In  an  early  study  by  Nychas,  Hershey 
and  Brodkey  (1973)  It  was  suggested  that  (what  appeared  to  be)  large 
transverse  convecting  vortices  in  the  outer  layer  were  somehow 
associated  with  the  bursting  events  near  the  wall.  Bakewell  and 
Lumley  (1967)  and  later  Blackwelder  (1978)  and  Blackwelder  and 
Eckel mann  (1979)  suggested  that  long  counter-rotating  streamwise 
vortex  pairs  are  an  important  feature  in  the  time-dependent  flow 
near  the  wall.  Falco  (1981,1982)  has  proposed  that  convected  ring¬ 
like  vortices  are  an  important  dynamical  feature  of  turbulent 
boundary  layers.  Perhaps  the  most  popular  vortex  model  is  the  hair¬ 
pin  vortex  which  has  been  proposed  as  a  basic  feature  of  wall -layer 
turbulence  by  Theodorsen  (1952),  Wlllmarth  and  Tu  (1967),  Kline 
et  al  (1967),  Offen  and  Kline  (1973,1975),  Smith  (1978),  Wallace 
(1982),  Head  and  Bandyopadhyay  (1981),  Perry,  Lim  and  Teh  (1981), 
Perry  and  Chong  (1982),  Smith  and  Metzler  (1983)  and  Acarlar  and 
Smith  (1984,  1987a,  1987b)  among  others. 

In  the  study  by  Nychas,  Hershey  and  Brodkey  (1973),  small 
solid  particles  were  suspended  in  a  water  flow  and  the  particle 
motions  were  photographed  using  a  high-speed  motion-picture  camera 
which  moved  in  the  flow  direction.  It  is  difficult  to  get  a  clear 
three-dimensional  picture  of  the  sequence  of  events  with  this 
approach.  Nevertheless,  Nychas  et  al  (1973)  observed  what  appeared 
to  be  transverse  vortices  convecting  in  the  outer  layer.  Such 
vortices  induced  a  decelerated  flow  near  the  wall  followed  by  a 
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region  of  acceleration.  It  was  concluded,  on  the  basis  of  the 
visual  observations,  that  ejection  events  at  the  wall  were  closely 
associated  with  the  presence  of  transverse  vortices  in  the  outer 
flow;  the  clear  suggestion  of  this  work  is  that  the  transverse 
vortices  induce  conditions  in  the  wall  region  that  cause  an 
ejection.  A  similar  mechanism  is  suggested  by  the  studies  of  Falco 
(1981,1982)  and  Falco  and  Wiggert  (1980);  in  these  experiments,  the 
flow  structure  of  the  turbulent  boundary  layer  was  visualized  by 
using  smoke  in  air  and  compact  vortex  structures,  which  appeared  to 
be  similar  to  vortex  rings,  were  observed;  the  ring-like  vortices 
were  believed  to  induce  eruptions  of  the  wall -layer  flow.  Detailed 
experiments  involving  vortex  rings  moving  toward  a  wall  by  Cerra  and 
Smith  (1982),  Walker,  Smith,  Cerra  and  Doligalski  (1987)  and  Didden 
and  Ho  (1986)  confirm  that  vortex  rings  moving  toward  a  well  do 
induce  eruptive  behavior  and  the  production  of  secondary  vortices 
through  an  unsteady  separation  effect  with  the  boundary  layer  near 
the  wall.  This  mechanism  will  be  discussed  in  more  detail  in 
Section  1.3. 

Another  picture  of  the  flow  near  the  wall  in  turbulent 
boundary  layer  has  been  proposed  by  Bakewell  and  Lumley  (1967); 
these  authors  examined  the  near-wall  region  of  a  pipe  flow  and  on 
the  basis  of  measurements  with  a  hot-film  anemometer  concluded  that 
the  most  energetic  velocity  fluctuations  were  associated  with  pairs 
of  counter-rotating  vortices.  Bakewell  and  Lumley  (1967)  suggested 
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that  the  dominant  structures  in  the  wall -layer  region  consist  of 
randomly  distributed  counter-rotating  vortices  which  are  elongated 
in  the  streamwise  direction.  A  similar  model  was  proposed  by 
Blackwelder  and  Eckelmann  (1979)  who  suggested  that  long  streamwise 
vortices  are  present  in  the  wall-layer  flow;  in  this  model,  the  rol 
of  the  counter-rotating  vortex  pairs  is  to  pump  low  momentum  fluid 
away  from  the  wall  and  to  lead  to  the  formation  of  the  low-speed 
streaks.  It  was  conjectured  that  the  process  terminates  with  the 
arrival  of  a  sweep  of  high  speed  flow  from  the  outer  region  which 
was  supposed  to  lead  to  a  highly  inflectional  streamwise  profile 
near  the  low-speed  streaks.  The  inflectional  velocity  profile  is 
believed  to  be  a  highly  unstable  situation  and  the  bursting  event 
was  considered  to  be  the  outcome  of  the  instability.  This  type  of 
model  is  purely  kinematical  and  assumes  a  great  deal  about  insta¬ 
bility  in  a  time-dependent  evolving  flow  field;  apart  from  these 
difficulties,  the  physical  origin  of  the  counter-rotating  vortex 
pairs  was  not  explained. 

A  third  type  of  vortex  structure  that  has  often  been  proposed 
as  a  basic  feature  in  turbulent  boundary-layer  flows  is  the  hairpin 
vortex.  A  variety  of  terms  have  been  used  to  describe  hairpin 
vortices;  Head  and  Bandyopadhyay  (1981)  have  pointed  out  that 
these  terminologies  really  imply  the  same  vortex  structure  but  at 
different  Reynolds  numbers.  The  classification  scheme  taken  from 
their  paper  is  reproduced  in  Figure  1.3, which  simply  implies  that 


Effect  of  Reynolds  number  on  eddies 

(a)  very  slow  Re  (vortex  loops); 

(b)  low-moderate  Re  (horseshoes); 

(c)  moderate-high  Re  (hairpins). 


at  higher  Reynolds  numbers  the  vorticity  tends  to  become  more  tightly 
concentrated  in  smaller  and  smaller  core  regions.  Head  and 
Bandyopadhyay  (1981)  carried  out  a  series  of  wind  tunnel  flow 
visualization  experiments  using  smoke;  their  results  indicate 
that  the  turbulent  boundary  layer  consists  of  a  "forest"  of  hairpin 
vortices  which  continually  stretch  in  shear  and  _ppear  to  be  inclined 
at  a  characteristic  angle  of  45°  to  the  wall.  A  similar  picture  has 
also  been  suggested  by  Perry  and  Chong  (1982)  who  conjectured  that 
the  structure  of  the  turbulent  boundary  layer  is  dominated  by  a 
systematic  heirarchy  of  hairpin  vortices.  There  are  some  appealing 
features  about  this  conjectured  type  of  structure  but  the  physical 
picture  appears  Incomplete.  For  example,  there  is  no  explanation 
concerning  the  processes  leading  to  the  production  of  new  vorticity 
from  the  wall  region;  In  addition  the  notion  of  a  series  of  hairpin 
vortices  which  continually  stretch  and  intertwine  in  shear  does  not 
appear  to  lead  to  a  complete  and  an  entirely  satisfactory  pattern  of 
evolution. 

Recently,  Smith  (1984)  has  proposed  a  conceptual  model  of  the 
physical  processes  that  occur  near  the  wall  in  a  turbulent  boundary 
layer  which  was  based  upon  extensive  flow  visualization  experiments 
as  well  as  a  previous  model  by  Offen  and  Kline  (1975).  In  the  latter 
study,  it  was  suggested  that  lift-up  of  the  low  speed  streaks  (which 
appeared  to  be  a  precursor  of  the  bursting  process)  was  associated 
with  the  temporary  imposition  of  an  adverse  pressure  gradient  on  the 


wall  layer  by  a  convected  disturbance.  Smith  (1984)  also  adopts 
this  view  and  suggests  that  the  effect  of  the  local  adverse  pressure 
gradient  Is  to  Induce  an  inflectional  behavior  in  the  streamwise 
velocity  profile  near  the  interface  of  the  streak  and  the  high  speed 
outer  region  fluid.  The  instantaneous  inflectional  profile  is 
assumed  to  be  unstable  and  to  give  rise  to  a  series  of  roll-ups  into 
a  number  of  hairpin  vortices.  As  the  hairpin  vortices  are  convected 
away  from  the  surface,  a  strong  vortex  stretching  process  in  the 
streamwise  direction  cause  intensification  in  the  legs  of  the  hair¬ 
pin  vortices;  thus, the  counter-rotating  legs  of  the  hairpin  vortices 
act  to  reinforce  the  streak,  thus  creating  an  appearance  of  streak 
persistence.  The  bursting  process  is  suggested  to  be  due  to  the 
breakup  of  the  low-speed  streak  which  is  lifted  by  the  counter¬ 
rotating  legs  of  the  hairpin;  rapid  ejection  may  be  due  to  the 
viscous  effects  induced  by  the  local  streamwise  adverse  pressure 
associated  with  the  moving  head  of  the  hairpin  vortex. 

In  order  to  observe  the  motion  of  hairpin  vortices  in  a  control 
led  environment,  Acarlar  and  Smith  (1984,  1937a,  1987b)  carried  out 
series  of  experiments  where  hairpin  vortices  were  artifically  gener¬ 
ated  in  a  laminar  boundary-layer  flow  on  a  plane  surface.  The  hair¬ 
pin  vortices  were  synthetically  generated  by  shedding  from  a  hemi¬ 
spherical  protuberance  on  the  plane  wall.  Flow  visualization  was 
carried  out  using  hydrogen  bubble  wires  and  it  was  determined  that 
convecting  hairpin  vortices  produce  flow  patterns  which  are  very 
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;  similar  to  those  observed  in  the  near-wall  region  of  a  turbulent 

> 

boundary  layer.  In  Figure  1.4,  which  is  reproduced  from  Acarlar  and 
Smith  (1984,  1987a),  the  bubble  line  patterns  generated  from  a  bub¬ 
ble  wire  placed  normal  to  the  wall  may  be  observed  for  (a)  a  convec- 

t 

j  ting  hairpin  vortex  in  a  laminar  boundary-layer  flow  and  (b)  in  a 

fully  turbulent  boundary  layer.  In  Figure  1.4,  the  quantity  X/R  = 

20  indicates  that  the  flow  pattern  is  at  a  streamwise  distance  equal 
I  to  20  radii  downstream  of  the  hemisphere;  Re.  is  the  Reynolds  number 

of  the  turbulent  boundary-layer  flow  based  on  momentum  thickness.  Ii 
Figure  1.5,  the  bubble  wire  was  placed  in  an  orientation  across  the 
span  of  the  flow  and  close  to  the  wall.  It  may  be  observed  that  the 
convecting  hairpin  vortices  produce  a  "streaky"  structure  in  the 
flow  near  the  wall  which  is  similar  to  the  typical  patterns  observed 
in  turbulent  boundary  layers.  In  addition  to  the  flow  visualization 
studies,  Acarlar  and  Smith  (1984,  1987a)  also  carried  out  quantita¬ 
tive  measurements  in  the  region  downstream  of  the  hemisphere;  it  was 
found  that  the  measurements  were  very  similar  to  those  typically 
obtained  in  a  fully  turbulent  boundary  layer. 

It  is  reasonable  to  conclude  from  the  experimental  evidence 
available  that  the  hairpin  vortex  is  a  potentially  important  flow 
structure  in  turbulent  flows  near  walls  and  worthy  of  further  study. 

1.3  Theoretical  Background 


1.3.1  Inviscid  Vortex  Dynamics 

In  recent  times,  there  has  been  an  increasing  interest  in  the 


(a)  Hairpin  vortex,  X/R  =  20 


(b)  Turbulent  boundary  layer.  Re  2200 

Figure  1.4  -  Side  view  conpari ‘.on  between  hairpin  vortex 
and  turbulent  boundary-1  aver  patterns. 


(a)  Hairpin  vortex,  X/R  =  20,  Yw-jrp/^  = 


development  of  numerical  algorithms  to  calculate  three-dimensional 
vortex  motion.  The  general  approach  may  be  described  as  follows. 

For  a  given  vortex  filament  and  a  given  initial  configuration,  a 
discretization  is  accomplished  by  placing  a  number  of  nodes  on  the 
filament;  some  approximation  to  the  Biot-Savart  integral  (Batchelor, 
1967)  is  then  used  to  compute  the  local  instantaneous  velocity 
induced  on  a  node  by  the  other  portions  of  the  vortex.  The  motion 
of  each  node  is  then  computed  numerically  thus  giving  the  position 
of  the  vortex  filament  at  the  next  time  step.  Such  calculations 
are  time  consuming,  notoriously  unstable  and  difficult  to  carry  out 
with  good  accuracy.  One  major  difficulty  is  that  a  vortex  filament 
having  a  core  of  zero  cross-sectional  area  induces  an  infinite 
velocity  on  itself  for  any  curved  portion  of  the  vortex;  the 
infinitesimally  thin  vortex  is  therefore  not  a  self-consistent 
dynamical  model  and  in  any  computational  algorithm  it  is  necessary 
to  take  into  account  the  fact  that  the  vortex  has  a  core  of  finite 
thickness.  In  addition  it  is,  in  principle,  necessary  to  account 
for  motions  within  the  core  and  how  such  motions  can  influence  the 
overall  motion  of  the  vortex. 

The  theoretical  problem  is  immediately  evident  when  the  Biot- 
Savart  integral  over  all  space  is  replaced  by  a  line  integral 
along  the  instantaneous  location  of  the  vortex  filament;  a  singular¬ 
ity  occurs  in  the  integrand  at  any  point  on  the  vortex  filament  where 
an  evaluation  of  the  velocity  induced  by  the  vortex  is  desired 


(this  location  will  be  subsequently  referred  to  as  the  field  point). 
This  singularity  arises  as  a  direct  consequence  of  neglecting  the 
thickness  of  the  vortex  core.  Over  the  years,  a  variety  of 
approaches  have  been  proposed  to  deal  with  this  difficulty  and 
these  will  be  discussed  in  more  detail  in  Section  2.3;  only  a  brief 
summary  will  be  given  here. 

The  first  attempts  at  the  calculation  of  vortex  motion  in 
three-dimensions  were  by  Hama  (1962)  and  Arms  and  Hama  (1965).  In 
their  "localized  induction  concept",  a  Taylor  series  was  used  to 
represent  the  integrand  near  the  field  point;  the  integration  was 
carried  out  only  over  a  short  portion  of  the  filament  and  the  con¬ 
tribution  of  the  rest  of  the  vortex  was  ignored.  Crow  (1970) 
introduced  the  "cut-off  method"  wherein  a  small  piece  of  the  Biot- 
Savart  integral  is  omitted  from  the  integration  thereby  bypassing 
the  singularity.  In  subsequent  studies,  Bliss  (1970)  and  Moore 
and  Saffman  (1972)  put  the  "cut-off"  procedure  on  a  firmer  basis 
by  selecting  the  cut-off  length  so  that  the  computational  algorithm 
agreed  with  known  exact  results  for  a  translating  circular  vortex 
ring.  The  "cut-off  method"  is  somewhat  inconvenient  from  a  computa¬ 
tional  point  of  view  and  an  alternative  (and  similar)  approach  was 
adopted  by  Moore  (1972).  In  this  technique,  a  small  parameter  is 
artiflcally  introduced  in  the  denominator  of  the  integrand  to  remove 
the  singularity;  the  small  parameter  is  selected  so  that  the  pro¬ 
cedure  gives  correct  results  for  the  circular  vortex  ring. 


Another  computational  approach  is  due  to  Leonard  (1975,1979) 
and  Widnall,  Bliss  and  Zalay  (1971).  In  this  technique,  the  self- 
induced  velocity  of  a  vortex  filament  is  computed  through  numerical 
integration  of  the  Biot-Savart  integral  except  at  the  field  point 
and  the  two  adjacent  nodal  points  on  the  vortex.  The  contribution 
due  to  these  three  points  is  estimated  by  fitting  a  vortex  ring 
locally  and  by  estimating  the  contribution  to  velocity  due  to  the 
arc  from  known  exact  results  for  the  vortex  ring.  Note  however  that 
Liu,  Krause  and  Ting  (1985)  point  out  that  this  process  does  not 
necessarily  give  the  correct  direction  of  self-induced  motion  of  the 
arc  surrounding  the  field  point. 

In  recent  times,  Callegari  and  Ting  (1978)  have  used  the 
method  of  matched  expansions  to  investigate  the  dynamics  of  vortex 
motion  in  three  dimensions  in  the  limit  of  vanishing  viscosity;  (see 
also  Bliss,  1970).  These  authors  were  able  to  determine  the  possible 
form  of  a  viscous  solution  in  the  vortex  core  and  to  match  this  to 
an  outer  inviscid  solution  described  by  the  solution  of  Biot-Savart 
integral.  It  was  demonstrated  that  the  apparent  singularity  in  the 
Biot-Savart  integral  is  removed  by  this  procedure  and  a  computational 
algorithm  was  developed  to  compute  the  evolution  of  vortex  motion. 
Note,  however,  that  the  method  at  present  is  restricted  to  closed 
vortex  loops. 


1.3.2  Influence  of  Convecting  Vortices  on  Viscous  Flow 
Near  Walls 


As  a  vortex  is  convected  over  a  wall,  it  induces  an  unsteady 
pressure  gradient  on  the  viscous  flow  near  the  wall.  The  effect  of 
the  vortex  motion  has  been  studied  by  Walker  (1978),  Doligalski 
et  al  (1980),  Doligalski  and  Walker  (1984),  Ersoy  and  Walker  (1985, 
1986)  and  Walker  et  al  (1987)  for  a  variety  of  two-  and  three- 
dimensional  boundary-layer  flows.  In  general,  these  studies  show 
that  a  moving  vortex  induces  a  region  of  adverse  pressure  gradient 
in  the  flow  near  th_*  wall  which  in  turn  drives  an  unsteady  boundary 
layer  separation  effect.  In  the  cited  studies,  it  emerges  that  a 
wide  variety  of  complex  separation  effects  occur;  however  all 
situations  have  one  feature  In  common.  Eventually  the  boundary- 
layer  flow  near  the  wall  develops  strong  outward  local  growth  which 
ultimately  leads  to  a  strong  inviscid-viscous  interaction  with  the 
outer  flow.  Specific  physical  situations  studied  are  as  follows. 

The  unsteady  boundary  layer  due  to  rectilinear  vortex  in  an 
otherwise  stagnant  fluid  above  a  plane  wall  was  considered  by 
Walker  (1978).  Inviscid  theory  predicts  that  the  vortex  will 
convect  at  constant  height  and  speed  in  the  velocity  field  of  the 
Image  vortex.  However,  a  short  time  after  the  initiation  of  the 
motion,  a  secondary  eddy  of  opposite  rotation  to  the  parent  vortex 
occurs  in  the  unsteady  boundary-layer  flow  near  the  wall.  The 
numerical  integrations  show  that  strong  and  accelerating  boundary- 
layer  growth  occurs  near  the  secondary  eddy  and  strongly  suggest 
that  the  boundary-layer  will  erupt  into  the  outer  inviscid  flow 


region.  The  experiments  of  Harvey  and  Perry  (1972)  confirm  that 
this  is  the  case  and  show  that  the  inviscid-viscous  interaction 
consists  of  the  ejection  of  the  secondary  eddy  from  the  boundary 
1 ayer . 

The  boundary  layer  due  to  a  two-dimensional  vortex  convected 
in  a  uniform  flow  and  a  shear  flow  above  a  plane  wall  was  considered 
by  Doligalski  and  Walker  (1984)  and  Doligalski  et  al  (1980) 
respectively.  At  low  convection  speeds,  the  formation  of  a  second¬ 
ary  eddy,  of  opposite  rotation  to  the  parent  vortex,  was  observed 
in  the  boundary- layer  flow  along  with  explosive  outward  boundary- 
layer  growth  near  the  eddy.  At  higher  convection  speeds,  secondary 
eddies  did  not  occur  but  a  narrow  band  of  strong  growth  was 
observed  in  the  calculations  which  was  conjectured  to  ultimately 
give  rise  to  an  interaction  with  the  outer  flow.  On  the  basis  of 
these  results  Doligalski  et  al  (1980)  proposed  a  regenerative 
mechanism  for  turbulent  flows  near  walls.  According  to  this  pro¬ 
cess,  the  viscous  flow  near  the  wall  begins  to  respond  to  the  motion 
of  a  convecting  parent  vortex  above  the  wall  giving  rise  to 
accelerating  growth  in  a  narrow  band  behind  the  moving  parent  vortex 
As  time  increases,  the  upwelling  fluid  penetrates  the  inviscid  flow 
region  above  the  wall  to  an  increasing  extent;  as  the  rising  fluid 
penetrates  the  cross  flow  above  the  wall,  it  is  ultimately  over¬ 
turned  in  a  strong  inviscid-viscous  interaction  leading  to  a  roll-up 
into  a  secondary  vortex. 


The  wall -layer  model  of  Blackwelder  and  Eckelmann  (1979)  is 
based  on  the  assumption  that  counter-rotating  vortex  pairs  play  an 
important  role  in  the  dynamics  of  the  near-wall  flow.  Ersoy  and 
Walker  (1985,1986)  have  computed  the  viscous  flow  development  near 
a  wall  due  to  counter-rotating  vortex  pairs  above  the  wall.  A  number 
of  cases  were  considered  and  although  a  variety  of  unusual  complex 
separation  phenomena  were  observed,  the  boundary- layer  flow  evolved 
toward  an  eruptive  state  in  all  cases.  When  the  sense  of  rotation 
of  the  vortices  is  such  that  the  flow  is  away  from  the  wall  between 
the  vortices,  a  boundary -layer  eruption  and  local  inviscid-viscous 
interaction  is  indicated  In  a  streamwise  region  between  the  parent 
vortex  centers;  when  the  sense  of  rotation  of  the  vortices  is 
reversed,  the  interaction  occurs  outboard  of  the  vortex  centers. 

One  feature  which  is  missing  from  the  aforementioned  studies 
and  which  is  important  in  the  three-dimensional  flows  observed  in 
turbulent  boundary  layers  is  vortex  stretching.  The  subsequent 
studies.  Walker  et  al  (1987)  and  Ersoy  and  Walker  (1985)  have 
examined  the  boundary  layer  induced  by  a  vortex  ring  moving  toward  a 
wall  and  a  vortex  loop  approaching  a  wall  at  angle  of  attack.  Again 
strong  boundary-layer  growth  was  observed  in  the  numerical  calcula¬ 
tions;  vortex  loops  stretch  as  they  approach  the  wall,  and  as  the 
local  vorticity  levels  increase,  the  eruptive  response  of  the 
boundary-layer  flow  appears  to  intensify. 
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1.4  Summary  and  Overview 

Current  experimental  research  suggests  that  the  motion  of 
vortices,  particularly  hairpin  vortices,  plays  an  important  part  in 
the  time-dependent  dynamics  of  the  turbulent  boundary  layer.  The 
studies  by  Acarlar  and  Smith  ( 1984, 1987a, 1987b)  suggest  that  hairpin 
vortices  are  regenerative  in  at  least  two  ways:  (1)  a  hairpin  in 
shear  appears  to  be  able  to  multiply  itself  to  the  side  in  the  sense 
that  secondary  hairpin  vortices  are  observed  to  form  in  the  spanwise 
direction  out  board  of  the  original  hairpin  vortex  and  (2)  the  crea¬ 
tion  of  secondary  hairpins  through  an  inviscid-viscous  interaction 
with  the  viscous  flow  near  the  wall  was  continually  observed.  In 
Figures  1.6(a)  and  1.6(b)  (Taken  from  Acarlar  and  Smith  (1984))  both 
of  these  processes  are  depicted  schematically.  In  the  present  study 
the  dynamics  of  hairpin  vortices,  with  particular  reference  to  the 
phenomena  observed  by  Acarlar  and  Smith  (1984),  will  be  investigated 
on  a  theoretical  basis. 

The  plan  of  this  study  is  as  follows.  In  Chapter  2,  a 
numerical  method  for  computing  the  motion  of  a  vortex  filament  which 
is  based  on  a  numerical  approximation  to  the  Biot-Savart  integral 
will  be  given.  It  was  determined  that  current  numerical  schemes  for 
vortices  having  small  core  radii  are  highly  unstable;  previous 
numerical  investigation  of  vortex  motion  (c.f.  Dhanak  and 
De  Bernardinis,  1981)  have  been  carried  out  for  an  unboundeu  fluid 
and  have  used  relatively  thick  core  sizes.  In  this  type  of  situa- 


tion  the  numerical  stability  problem  is  not  evident.  In  the  present 
study,  it  was  important  to  develop  a  reliable  calculation  method  for 
vortices  with  small  cores;  this  is  because  the  main  interest  was  in 
flows  with  walls  and  once  the  vortex  core  touches  the  wall,  the 
Biot-Savart  calculation  must  terminate. 

In  Chapter  3,  the  numerical  algorithm  is  tested  by  computation 
of  the  trajectory  of  closed  vortex  loops  near  a  wall.  In  addition, 
a  series  of  calculations  is  performed  to  compute  the  trajectories  of 
simulated  "hydrogen  bubble  lines"  in  the  flow  as  the  vortex  passes 
by.  These  computations  show  that  a  moving  vortex  loop  causes  the 
simulated  "bubble  lines"  to  form  into  a  streaky  type  of  structure 
near  the  wall.  The  indication  of  these  results  is  that  a  convected 
vortex  loop  is  a  potential  creator  of  wall -layer  streaks  in  a 
turbulent  boundary  layer. 

In  Chapter  4,  numerical  solutions  are  carried  out  to  show  the 
evolution  of  a  hairpin  vortex  in  both  a  stagnant  flow  and  a  shear 
flow.  The  evolution  of  the  vortex  in  a  shear  flow  is  consistent  with 
the  experimental  results  of  Acarlar  and  Smith  (1984).  Simulations  of 
hydrogen  bubble  lines  also  show  that  the  hairpin  vortex  is  a  poten¬ 
tial  creator  of  wall-layer  streaks.  Finally,  a  numerical  calculation 
is  carried  out  to  show  the  interaction  that  occurs  between  two  hair¬ 
pin  vortices  convecting  in  a  shear  flow. 

In  Chapter  5,  a  method  is  described  which  allows  the  determina¬ 
tion  of  the  viscous  response  of  the  boundary-layer  flow  near  the  wall 


to  the  motion  of  the  convected  hairpin.  The  unsteady  boundary- 
layer  flow  is  three-dimensional  and  unsteady  and  the  computation  of 
this  flow  would  require  computer  resources  well  beyond  what  was 
available  in  the  course  of  this  study.  Consequently  a  method  was 
devised  which  allows  the  calculation  of  the  boundary-layer  evolution 
along  the  symmetry  plane  of  the  hairpin  vortex.  The  computed 
results  are  strongly  suggestive  that  the  moving  hairpin  vortex 
rapidly  evokes  an  eruptive  response  in  the  viscous  flow  near  the 
wall.  Detailed  conclusions  are  discussed  in  Chapter  6. 
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CHAPTER  2 

THE  VORTEX  CALCULATION  METHOD 

2.1  Introduction 

Three-dimensional  vortex  motions  play  an  important  role  in  a 
variety  of  practical  situations  involving  fluid  flow  and  consequent¬ 
ly  several  computational  methods  have  been  developed  over  the  past 
two  decades  to  predict  the  evolution  of  three-dimensional  vortex 
motion.  In  situations  where  the  vorticity  in  the  flow  field  is  con¬ 
centrated  in  thin  vortex  filaments  which  convect  in  an  otherwise 
irrotational  flow,  calculation  procedures  are,  in  general,  based  on 
the  Biot-Savart  law  and  are  Langrangian  In  nature;  the  vortex  fila¬ 
ment  is  discretized  at  some  Initial  time  using  a  number  of  nodal 
points  on  the  space  curve  defining  the  filament  and  the  subsequent 
motion  of  the  nodal  points  is  evaluated  as  the  filament  distorts. 

In  this  chapter,  a  calculation  method  for  the  evolution  of  vortex 
filaments  Is  developed.  The  general  features  of  the  Biot-Savart  law 
are  discussed  in  section  2.2  and  some  early  calculation  schemes  are 
reviewed  in  section  2.3.  In  this  study  an  adaptation  of  a  procedure 
due  to  Moore  (1972)  was  used  to  compute  the  evolution  of  vortices 
and  in  section  2.4,  a  basic  numerical  Instability  in  this  scheme  is 
described.  An  improved  calculation  procedure  Is  developed  in 
section  2.5. 

2.2  The  Biot-Savart  Law 

Consider  a  Cartesian  coordinate  system  (x,y,z)  with  origin  at 
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0  and  let  the  velocity  vector  at  a  point  p0(x0>y0»zQ)  be  denoted  by 
u(XQ)  where  XQ  Is  the  vector  from  the  origin  to  the  point  P  .  In 
general  (Batchelor,  1967),  the  velocity  at  PQ  is  related  to  the 
vorticlty  distribution  Z  throughout  space  by  the  Biot-Savart  law 
which  Is  given  by. 


(X0-X)  X  >ji 

|x0-xp 


dv(x)  4  sext(x0) 


(2.1) 


In  equation  (2.1)  the  volume  Integral  is  carried  out  over  all  space 
where  the  vorticity  is  non-zero  and  uext  is  an  background  external 
Irrotational  velocity  field.  The  variable  t  is  an  Integration  vari- 
able  and  the  vector  X  -  XQ  is  directed  from  the  field  point  PQ  to 
locations  in  the  volume  V  as  the  integration  in  equation  (2.1)  is 
carried  out. 


Now  consider  the  physical  situation  where  the  majority  of  the 
flow  field  is  Irrotational  and  the  vorticity  is  concentrated  in 
vortex  tubes  of  small  cross-section.  For  a  single  vortex  tube,  as 
the  cross-sectional  area  of  the  tube  approaches  zero,  the  vortex 
approaches  being  a  line  filament  defined  by  a  space  curve.  The 
vorticity  vector  at  any  point  on  the  curve  may  be  written  as. 


Z  =  nt  ,  (2.2) 

A 

where  t  is  the  unit  tangent  to  the  curve  (as  indicated  in  Figure  2.1) 
and  n  is  the  local  magnitude  of  the  vorticity.  Note  that  the  direc- 
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Figure  2.1  -  A  portion  of  a  vortex  filament 
(Shown  for  r<0). 
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tlon  of  positive  vorticity  is  determined  by  the  right  hand  rule  and 
for  the  situation  in  Figure  2.1,  a  vortex  having  the  indicated  sense 
of  rotation  has  negative  n.  Since  the  vorticity  is  assumed  to  be 
non-zero  only  within  the  small  vortex  core,  the  integrand  in  equa¬ 
tion  (2.2)  vanishes  except  within  the  core  where  the  volume  element 
may  be  written. 


dV  =  dAcd£  .  (2.3) 

Here  ds.  and  dAc  are  differential  length  and  cross-sectional  area 
elements  along  the  core,  respectively.  The  integral  of  the  vorticity 
distribution  over  the  cross-sectional  area  of  the  core  may  be  related 
to  the  circulation  about  the  core  by  Stokes'  theorem  and  it  is  easily 
shown  that 
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jf  ttdAc  =  ^  u  •  ds 


r  . 


(2.4) 


Here  C'  Is  any  closed  contour  surrounding  the  vortex  core  and  ds  is  a 
length  element  on  this  path.  According  to  Kelvin's  theorem 
(Batchelor,  1967),  the  circulation  r  is  constant  in  time  and  also 
constant  along  the  length  of  the  vortex  filament.  Consequently,  the 
Biot-Savart  law  (2.1)  becomes. 


i(x0)=-^/^e 

0  4*c  X.-X  3 


x  dX  *  Vt(X0)  . 


(2.5) 
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where  dX  is  a  differential  length  element  along  the  space  curve  C 

-y  -y 

defining  the  current  position  of  the  vortex.  Since  u(XQ)  =  3XQ/9t, 
equation  (2.5)  apparently  defines  a  differential  equation  from 
which  the  evolution  of  each  point  along  the  vortex  can  be  traced 
forward  in  time.  The  difficulty  is  that  for  XQ  on  the  curve  C,  the 
integrand  is  singular  and  consequently  the  velocity  field  described 
by  equation  (2.5)  must  be  interpreted  as  an  outer  inviscid  solution 
which  is  not  uniformly  valid  near  the  core  of  the  vortex.  To  de¬ 
fine  a  rational  calculation  procedure,  it  is  in  theory  necessary  to 
match  any  solution  obtained  from  equation  (2.5)  with  an  inner  solu¬ 
tion  in  the  vortex  core  and  thereby  take  into  account  how  the  evolv¬ 
ing  flow  in  the  vortex  core  influences  the  overall  motion  of  the 
vortex.  It  emerges  that  this  is  a  rather  complex  problem  and 
although  recent  progress  has  been  made  (Callegari  and  Ting  (1978))  in 
this  area,  a  variety  of  approximate  methods  have  been  used  over  the 
years.  Some  of  these  techniques  are  discussed  in  detail  by  Ersoy 
and  Walker  (1985b)  and  only  a  brief  summary  is  given  here. 

2.3  Previous  Calculation  Methods 

For  vortex  motions  in  three-dimensions,  there  is  one  flow 
situation  where  the  core  solution  can  be  determined  and  subsequently 
matched  to  the  outer  inviscid  flow  field  described  by  the  Biot- 
Savart  integral;  this  Is  the  problem  of  a  circular  vortex  ring  which 
propogates  into  an  otherwise  Irrotational  flow  and  has  been  con¬ 
sidered  by  Tung  and  Ting  (1967)  and  Saffman  (1970).  Consequently 


there  is  one  exact  solution  for  three-dimensional  flows  and  this 
solution  has  been  used  to  structure  a  number  of  vortex  calculation 
methods.  In  the  early  calculation  methods  of  Widnall ,  Bliss  and 
Zalay  (1971)  and  Leonard  (1975,1979)  the  curve  defining  the  vortex 
filament  was  first  discretized  into  a  number  of  nodal  points;  as 
the  vortex  filament  moves,  the  velocity  induced  by  the  filament  at  a 
particular  nodal  point  was  calculated  by  a  numerical  integration  of 
the  Biot-Savart  law  except  over  the  segments  of  the  vortex  on  either 
side  of  the  nodal  point  in  question.  The  contribution  of  the  in¬ 
duced  velocity  due  to  the  segments  of  the  vortex  adjacent  to  the 
nodal  point  was  then  estimated  by  fitting  a  vortex  ring  through  the 
three  adjacent  points;  the  velocity  Induced  by  this  arc  at  the  nodal 
point  Is  then  estimated  using  the  known  results  for  a  circular 
vortex  ring.  The  approach  has  been  termed  a  "patched"  solution 
technique  by  Liu,  Tavantzis  and  Ting  (1984)  and  does  not  neces¬ 
sarily  give  the  correct  direction  of  the  self-induced  motion  locally. 

Another  approach  to  the  cut-off  method  originally  described  by 
Crow  (1970)  in  his  Investigation  of  the  stability  of  aircraft  trail¬ 
ing  vortices.  In  this  method,  a  portion  of  the  filament  near  the 
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singularity  at  XQ  is  simply  omitted  In  the  Biot-Savart  integral. 

Let  s  measure  distance  along  the  space  curve  defining  the  instantan¬ 
eous  vortex  position  and  let  sQ  denote  the  value  of  s  corresponding 

~y 

to  XQ;  in  the  cut  off  method,  the  approximation  to  the  Biot-Savart 
integral  is. 
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The  notation  in  equation  (2.6)  implies  that  an  interval  of  length  6c 
has  been  omitted  in  the  integration  on  both  sides  of  s  =  sQ;  Crow 
(1970)  determined  a  value  for  the  cut-off  length  6  by  comparing 
the  result  predicted  by  equation  (2.6)  when  the  vortex  is  a  circular 
ring  with  the  known  exact  solution.  As  a  result  of  this  comparison, 
a  value  of  &c  =  0.642  a  was  suggested  by  Crow  (1970);  here  a  denotes 
the  small  (not  nonzero)  radius  of  the  vortex  core.  Moore  and 
Saffman  (1972)  were  subsequently  able  to  give  a  formula  for  6c 
according  to. 


log{-^-}  =  |  -  ^/v'2rdr  +  77- J  w,2rdr 
0  0 


(2.7) 


where  r  measures  distance  from  the  center  of  the  vortex  core  and  v1 
and  w'  denote  the  swirl  and  axial  velocities  In  the  core.  The  right 
side  of  equation  (2.7)  was  obtained  from  detailed  results  for  the 
flow  in  the  core  of  a  circular  vortex  ring;  note  that  the  result  in 
equation  (2.7)  relates  the  cut-off  length  to  quite  general  velocity 
distributions  within  the  core  but  also  that  the  validity  of  the  cut 
off  method  is  assumed. 

From  a  computational  standpoint,  the  cut-off  method  of  Crow 
(1970)  is  normally  not  a  convenient  algorithm  for  the  following 
reason.  Calculations  involving  the  Biot-Savart  law  are  Lagrangian 


w 


in  character;  the  space  curve  defining  the  vortex  at  some  initial 
instant  is  discretized  into  a  number  of  nodal  points  which  then  may 
be  considered  to  be  Lagrangian  markers.  The  object  of  the  numerical 
integration  is  to  trace  the  path  of  these  markers  and  thus  determine 
the  shape  of  the  vortex  filament  at  subsequent  times.  To  use  Crow's 
(1970)  method,  it  is  necessary  to  subtract  a  portion  of  the  integra¬ 
tion  range  near  each  point  sQ  in  equation  (2.6);  generally  the  cut¬ 
off  length  will  not  correspond  to  the  mesh  length  and  indeed  will 
change  with  time  in  general.  An  alternative  approach  is  due  to 
Moore  (1972)  and  has  its  origins  In  a  study  by  Rosenhead  (1930);  in 
this  technique^the  singularity  in  the  integrand  is  removed  by 
artificially  inserting  a  small  parameter  u  into  the  denominator  and 
equation  (2.5)  is  written  according  to, 


u(sQ,t)  = 
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2+u2} 


372 


ds  + 


uext(so> 


(2.8) 


Here  s  Is  a  coordinate  along  the  vortex  and  sQ  denotes  a  specific 
location  on  the  vortex.  The  integral  is  now  taken  over  the  entire 
vortex  and  an  expression  for  the  small  parameter  u  is  obtained  by 
Insisting  that  equation  (2.8)  give  correct  results  for  the  circular 
vortex  ring.  The  following  expression  is  obtained  for  u  (Moore, 
1972;  Ersoy  and  Walker,  1985b). 


v '  ?r dr 


+  Y?-  /  w':rdr 


(2.9) 
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It  may  be  observed  from  equations  (2.8)  and  (2.9)  that  the  internal 
structure  of  the  vortex  core  flow  influences  the  overall  motion  of 
the  vortex  only  through  the  parameter  p  In  equation  (2.8)  which  is 
calculated  as  an  integrated  effect  over  the  core  via  equation  (2.9). 

It  might  be  expected  a  priori  that  the  core  radius  a  could  be  a 

function  of  position  s  along  the  vortex  as  well  as  of  t.  However, 
Moore  and  Saffman  (1972)  have  argued  that  the  vortex  core  does  not 
vary  in  the  axial  distance  and  is  uniform  at  any  instant  in  time; 
the  physical  reasoning  is  that  any  local  non-uniformities  in  the 
vortex  core  will  be  rapidly  smoothed  out  by  travelling  internal 
waves  on  a  time  scale  which  is  much  smaller  than  the  time  scale 

associated  with  the  overall  motion  of  the  vortex  itself.  Thus  the 

vortex  core  radius  may  be  regarded  as  a  function  of  time  alone  and 
Moore  and  Saffman  (1972)  also  show  that  , 

L(t)a2(t)  =  constant  ,  (2.10) 


where  L(t)  is  the  total  arc  length  of  the  closed  vortex  filament  at 
time  t;  this  result  simply  states  conservation  of  volume  of  the 
vortex  core.  It  follows  from  equation  (2.10)  that  the  core  radius  at 
any  instant  in  the  motion  is  described  by, 


a(t)  =  a  {j-  i 
o  c 


ii 

as 


(2.11) 
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where  aQ  and  Lq  are  the  initial  core  radius  and  vortex  length, 
respectively.  It  is  worthwhile  to  note  that  these  are  rather  strong 
results  which  have  been  confirmed  by  the  asymptotic  analysis  of 
Callegari  and  Ting  (1978).  However  the  results  apply  only  to  a 
vortex  filament  which  is  closed  and  therefore  of  finite  length.  In 
addition,  it  appears  that  the  theoretical  description  is  incomplete 
since  it  does  not  account  for  situations  where  local  variations  in 
core  size  are  observed  before  a  phenomenon  known  as  vortex  break¬ 
down  occurs. 

Moore  (1972)  has  used  equation  (2.8)  to  compute  the  motion  of 
aircraft  trailing  vortices;  the  method  was  also  subsequently  used  by 
Dhanak  and  De  Bernardinis  (1981)  and  by  Dhanak  (1981)  to  calculate 
the  evolution  of  elliptic  vortex  rings  and  a  vortex  filament 
approaching  a  rigid  sphere, respectively.  In  these  calculations, 
typical  values  of  the  vortex  core  radius  used  were  on  the  order  of 
0.2.  Although  the  integrand  in  equation  (2.8)  is  finite  everywhere, 
it  may  become  quite  large  near  sQ,  particularly  in  situations  where 
the  vortex  is  undergoing  stretching  and  a  (and  hence  u)  decreases.  To 


avoid  potential  inaccuracies  in  the  numerical  evaluation  of  the 


integral  near  sQ,  Moore  (1972)  obtains  the  leading  term  in  an 


expansion  around  s  =  sQ  according  to, 


3X 


(Vx) 


3s  X  { lx  -XI ?+u?}3/2 


#0  * 


■  32X: 


(2.12) 


*\i 


«s-so>  <5iV“2> 


and  the  subscript  o  indicates  that  the  partial  derivatives  are 
evaluated  at  s  =  sQ;  the  leading  order  form  of  the  integrand  near 
s  =  sQ  is  added  and  subtracted  to  the  integrand  to  obtain. 


u(s0,t) 
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The  apparent  advantage  of  this  procedure  is  that  the  integrand  in 
the  first  integral  is  0(1)  everywhere  while  the  second  integral  may 
be  evaluated  analytically.  The  method  was  used  with  good  success  by 
Moore  (1972),  Dhanak  and  De  Bernardinis  (1981)  and  Dhanak  (1981); 
however,  in  all  of  these  studies  rather  large  values  of  the  core 
radii  (a-0.1-0.3)  were  used.  In  the  present  study,  one  objective 
was  to  compute  the  evolution  of  hairpin  vortices  close  to  solid 
walls  and  to  avoid  situations  where  the  vortex  core  impacted  the 
wall  at  relatively  early  stages,  It  was  considered  desirable  to  be 
able  to  compute  cases  with  very  small  vortex  cores.  It  emerged 
that  the  method  of  Moore  (1972)  exhibited  strong  numerical  insta¬ 
bility  for  small  values  of  a;  a  modified  formula  based  on  equation 
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(2.8)  which  is  suitable  for  small  a  was  developed  in  this  study 
and  this  will  be  discussed  in  section  2.5. 


Another  approach  to  the  computation  of  vortex  motion  in  three- 
dimensions  was  originally  described  by  Arms  and  Hama  (1965)  and  this 
is  the  "localized  induction  approximation".  In  this  approach,  it 
is  assumed  that  the  dominant  contribution  to  the  self-induced 
vortex  velocity  occurs  locally;  in  this  approach,  the  integrand  in 
the  Biot-Savart  law  is  approximated  everywhere  by  its  asymptotic 
form  near  s  =  sQ  and  it  may  be  shown  that, 


(M)  x  (i£X) 

u(s  )  -  -iL9-_.-Ls._o  1  (1}  +0(1)  .  (2.15) 

|/lXx  3  a 

| ' 3S'0 

Note  that  this  result  may  be  obtained  from  equation  (2.14)  by 
ignoring  the  first  integral  and  integrating  P(s)  from  s  =  -i  to  +z 
where  z  is  an  0(1)  number;  the  integral  over  the  image  as  well  as 
the  external  field  are  also  ignored.  In  equation  (2.15),  a  is 
assumed  to  be  small  so  that  -1 og ( a ) > >1  and  the  first  term  of  equa¬ 
tion  (2.15)  is  presumed  to  describe  the  leading  order  behavior  of 
the  velocity  field  everywhere. 

The  "localized  induction  approximation"  was  used  by  Arms  ai.d 
Hama  (1965)  and  more  recently  by  Aref  and  Flinchem  (1984)  to  obtain 
some  interesting  results.  In  the  latter  study,  the  term  log(l/a) 
was  regarded  as  a  large  but  unknown  quantity  which  could  be  absorbed 
in  a  time  scale;  results  which  bear  a  similarity  to  those  obtained 


in  the  present  study  were  found  for  the  propogation  of  a  three- 
dimensional  distortion  on  a  two-dimensional  vortex.  However,  it  was 
evident  at  the  outset  that  three-dimensional  vortex  motion  involves 
some  rather  complex  phenomena  which  often  cannot  be  obtained 
through  physical  intuition;  in  addition,  the  effect  of  the  wall  was 
of  prime  r.^cern.  For  these  reasons,  it  was  considered  important 
to  carry  out  numerical  integrations  of  the  full  Biot-Savart  law  and 
to  account  for  the  presence  of  the  image  vortex. 

One  other  approach  should  be  noted.  Recently  Liu,  Tavantzis 
and  Ting  (1984)  have  described  a  computational  method  for  closed 
vortex  loops  which  is  potentially  more  general  than  some  of  the 
methods  described  here.  The  computational  method  is  based  on  a 
previous  analytical  study  by  Callegari  and  Ting  (1978)  in  which 
asymptotic  methods  were  used  to  write  the  Biot-Savart  integral  in 
equation  (2.1)  as  a  part  which  is  bounded  everywhere  (the  finite 
part)  plus  a  part  which  becomes  singular  as  the  curve  defining  the 
vortex  core  is  approached;  the  analysis  is  valid  for  small  vortex 
cores.  Within  the  core,  the  flow  is  viscous  and  Callegari  and  Ting 
(1978)  investigate  possible  unsteady  solutions  of  the  core  equa¬ 
tions;  these  are  matched  to  the  outer  inviscid  solution  corresponding 
to  the  singular  part  of  the  Biot-Savart  integral.  Consequently,  in 
this  approach  the  evolving  flow  in  the  vortex  core  influences  the 
overall  motion  of  the  vortex  loop  to  an  increasing  extent  as  time 
passes.  The  details  of  the  asymptotic  analysis  are  complex  and 
will  not  be  summarized  here.  To  date,  the  method  has  not  been 
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extended  to  three-dimensional  vortex  filaments  of  infinite  length, 
which  are  of  principal  interest  in  this  study. 

2.4  An  Instability  Associated  With  the  Moore-Rosenhead  Method 
In  the  present  study,  the  calculations  of  three-dimensional 
vortex  motion  were  based  on  the  Moore-Rosenhead  model  equations  (2.8) 
and  (2.9)  for  the  Biot-Savart  integral.  The  algorithm  based  on  this 
model  is  convenient  for  computational  purposes  and  is  believed  to 
realistically  model  the  evolution  of  vortex  motion.  In  order  to 
avoid  loss  of  accuracy  in  evaluation  of  the  integrand,  the  procedure 
suggested  by  Moore  (1972)  given  in  equation  (2.14)  was  used.  How¬ 
ever  after  running  several  test  cases,  it  was  determined  that  the 
scheme  suffered  from  serious  instabilities  which  generally  became 
worse  with  either  decreasing  core  size  or  an  increased  number  of 
points  on  the  vortex.  To  illustrate  the  problem,  consider  the  case 
of  a  circular  vortex  ring  which  is  moving  on  a  normal  trajectory 
toward  a  plane  wall;  the  exact  solution  for  the  problem  is  known 
(Walker,  Smith,  Cerra  and  Do! i gal  ski ,  1987)  and  it  provides  a  con¬ 
venient  test  case  for  the  present  calculation  methods.  Since, 

■  (2-,6) 

equation  (2.14)  defines  an  initial  value  problem;  for  a  specified 
number  of  mesh  points  along  the  vortex,  equation  (2.14)  may  be 
integrated  forward  in  time  to  track  the  motion  of  each  individual 


point  (denoted  by  sQ)  on  the  vortex.  The  results  of  such  a  calcula¬ 
tion  are  shown  in  Figure  2.2;  here  a  vortex  ring  (of  radius  1)  has 
been  started  at  a  height  0.5  units  away  from  the  wall.  The  vortex 
core  was  assumed  to  be  in  solid  body  rotation  with  core  radius 
a  =  0.005  and  with  v1  =  rr/2TTa2,  w*  =  0  it  follows  from  equation 
(2.9)  that, 


In  the  calculation  depicted  in  Figure  2.2,  there  were  100  nodal 
points  on  the  vortex  and  the  time  increment  was  At  =  0.005.  It  may 
be  observed  that  the  ring  proceeds  smoothly  toward  the  wall  but  that 
wiggles  have  developed  in  the  ring  by  t  =  0.05.  Once  "wiggles" 
develop  in  the  ring  shape  the  calcuation  rapidly  breaks  down  and 
gives  nonsensical  results  as  indicated  at  t  =  0.055.  The  onset  of 
the  appearance  of  the  wiggles  can  be  delayed  somewhat  by  reducing 
the  time  step  but  the  problem  is  made  worse  by  using  more  mesh  points 
to  discretize  the  vortex. 

The  difficulty  depicted  in  Figure  2.2  is  due  to  a  deficiency  in 
Moore's  (1972)  method  and  has  not  been  previously  observed  by  Moore 
(1972),  Dhanak  and  De  Bernardinis  (1981)  and  Dhanak  (1981)  because 
these  authors  have  used  relatively  thick  cores  (a=0.1  to  0.35). 

Here,  because  of  the  presence  of  the  wall,  it  was  considered 
desirable  to  be  able  to  use  very  small  cores  and  consequently  it  was 
necessary  to  find  a  resolution  of  the  problem. 


The  difficulty  which  gives  rise  to  the  instability  is 
associated  with  the  fact  that  for  small  values  of  a  (and  hence  u)  the 
demoninator  in  the  first  integrand  in  equation  (2.14)  is  small  near 
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s  =  sQ  and  hence  inaccuracies  in  the  derivative  sX/as  are  magnified. 
The  integrand  in  the  first  integral  in  equation  (2.14)  is  plotted  in 
Figure  2.3(a)  for  a  typical  point  s  at  t  =  0.01  for  the  calculation 
depicted  in  Figure  2.2.  It  may  be  observed  that  the  integrand 
exhibits  anamolous  behavior  near  s  =  sQ;  it  does  not  tend  smoothly 
to  zero  and  shows  two  small  negative  spikes  on  either  side  of  s  =  sQ. 
If  the  number  of  mesh  points  is  increased  to  200,  the  size  of  the 
spikes  diminishes  but  the  phenomena  persists.  With  increasing  mesh 
points  on  the  vortex,  the  numerical  evaluation  of  the  gradients 
aX/as  is  in  general  more  accurate  and  this  is  reflected  in  the 
diminished  levels  of  the  spikes  in  Figures  2.3.  In  Appendix  A,  it 
is  verified  that  the  problem  is  due  to  a  inaccurate  evaluation  of 
the  first  integrand  in  equation  (2.14);  a  basic  modification  of  the 
scheme  is  considered  next. 

2.5  The  Present  Method 

The  motivation  behind  the  algorithm  (2.14)  described  by  Moore 
(1972)  was  to  obtain  an  accurate  evaluation  of  the  integrand  near 
s  =  sQ.  However  the  results  of  the  previous  section  and  the 
discussion  in  Appendix  A  shows  that  there  are  still  difficulties 
with  the  procedure  for  vortices  with  small  cores.  To  eliminate 
these  problems,  the  integrand  in  the  first  integral  of  equation 


(2.14)  was  represented  by  higher  order  terms  in  an  expansion  about 


s  =  sQ;  the  details  are  described  in  Appendix  B  and  if  R(s,sQ) 
denotes  the  integrand  in  equation  (2.14)  according  to, 
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where  P(s)  is  defined  by  equation  (2.13),  it  can  be  shown  that 


R(s,s0)  =  D3(s-s0)3  +  D4(s-sq)4  +  — -  ,  (2.19) 


where  and  are  vectors  defined  by  equations  ( B .9 )  and  (B.10)  in 
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Appendix  B  in  terms  of  vector  products  of  derivatives  of  X  evaluated 
at  s  =  s0. 

In  Figure  2.4,  the  exact  value  of  the  integrand  R(s,sQ)  is  plot¬ 
ted;  here  the  expressions  for  a  ring  given  in  equations  (A. 2)  and 
(A. 3)  in  Appendix  A  have  been  used  to  evaluate  the  gradients  in 
equation  (2.18).  The  result  for  R(s,sQ)  using  numerical  differentia¬ 
tion  to  compute  the  gradients  and  the  local  approximation  (2.19)  is 
also  plotted  in  Figure  2.4.  It  may  be  observed  that  the  local 
approximation  (2.19) represents  R  very  well  for  small  (s-sQ)  and  that 
the  behavior  depicted  In  Figure  2.3  no  longer  occurs. 

A  calculation  method  for  the  vortex  ring  may  now  be  structured 
as  follows.  In  equation  (2.14),  the  first  Integral  is  split  into 
three  parts  and. 


Comparison  of  local  approximation  to  the 
integrand  R(s,s  )  to  the  exact  answer  for 
a  ring  vortex.  0 
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In  the  ranges,  [l,h),  [-1  ,-h) ,  the  first  integrand  is  evaluated  as 
before  while  in  the  interval  [-h,h]  surrounding  the  point  s  =  sQ, 
the  integrand  is  evaluated  using  the  series  representation  (2.19). 
Calculations  were  once  again  carried  out  for  the  vortex  ring  moving 
toward  the  wall  (r<0)  and  the  results  are  shown  in  Figure  2.5. 

Here  a  total  of  100  mesh  points  were  used  on  the  ring  with  an 
initial  vortex  core  radius  of  0.005  and  at  time  step  of  0.002.  It 
may  be  observed  that  the  calculation  proceedes  without  difficulty 
until  the  ring  is  quite  close  to  the  wall  in  contrast  to  the 
results  of  the  computation  depicted  in  Figure  2.2.  Note  that  in 
equation  (2.20),  the  value  of  h  must  be  selected  small  enough  so 
that  the  series  in  equation  (2.19)  is  convergent;  for  the  calcula¬ 
tions  depicted  in  Figure  2.5  h  =  3as  and  the  integrals  were  evalua¬ 
ted  by  Simpson's  rule. 


CHAPTER  3 


THE  INFLUENCE  OF  MOVING  VORTEX  LOOPS 
3.1  Introduction 

It  is  commonly  believed  that  vortex  motions  play  an  important 
role  in  the  dynamics  of  the  turbulent  boundary  layer.  Both  closed 
loop  vortices  (rings)  [Falco  (1982,1983),  Falco  and  Wiggert  (1980)] 
and  hairpin  vortices  [Acarlar  and  Smith  (1984, 1987a, 1987b)]  have  been 
proposed  as  principal  vortex  structures  in  the  outer  region  of  the 
turbulent  boundary  layer.  Wall- layer  streaks  are  observed  in  the 
near-wall  region  of  turbulent  flows  during  the  relatively  long 
periods  of  time  when  the  wall -layer  flow  is  not  strongly  interacting 
with  the  essentially  inviscid  flow  away  from  the  wall;  (Walker, 
Scharnhorst  and  Weigand,  1986).  It  is  therefore  of  interest  to 
understand  the  types  of  vortex  motion  which  will  induce  a  streaky 
structure  in  the  flow  near  a  wall.  In  this  chapter,  the  effect  of 
convecting  vortex  loops  will  be  considered;  the  question  of  con- 
vecting  hairpin  vortices  will  be  addressed  in  Chapter  4. 

The  numerical  calculations  that  will  be  described  here  may  be 
thought  of  as  an  inviscid  simulation  of  the  motion  of  hydrogen 
bubble  lines.  First  the  trajectory  of  a  vortex  ring  was  computed 
using  the  method  developed  in  Chapter  2;  for  situations  where  the 
ring  is  initially  inclined  at  an  angle  to  the  wall,  the  ring 
rapidly  distorts  into  a  vortex  loop  due  to  the  effect  of  the  wall. 

As  the  vortex  loop  convects  above  the  wall,  it  is  of  interest  to 


understand  the  nature  of  the  induced  flow  near  the  wall;  in 
typical  flow  visualization  experiments,  hydrogen  bubble  lines  are 
used  to  mark  the  flow  and  the  nature  of  the  flow  is  then  inferred 
from  the  subsequent  motion  of  the  bubble  lines.  In  the  present 
study,  the  bubble  lines  were  simulated  by  numerically  tracking  the 
motion  of  several  initially  straight  lines  near  the  wall  as  the 
vortex  loop  convected  past.  The  objective  of  this  simulation  is 
to  determine  whether  the  convecting  vortex  will  induce  the  marker 
lines  to  form  a  streaky  structure. 

3.2  The  Trajectory  of  Vortex  Loops 

Consider  a  vortex  ring  having  an  initial  radius  R  which  is 
initially  inclined  at  an  angle  a  to  the  wall  as  indicated  in 
Figure  3.1;  here  the  a  measures  the  angle  between  the  normal  to 
the  plane  of  the  ring  and  the  wall.  If  a  =  tt/2,  the  ring  will 
remain  circular  and  either  move  toward  or  recede  from  the  wall 
depending  on  the  sense  of  the  circulation;  as  the  ring  moves  to¬ 
ward  the  wall,  the  ring  radius  will  increase.  If  a  /  tt/2,  the 
ring  will  immediately  start  to  change  into  a  "rubber-band"  loop 
shape.  In  this  section,  the  results  for  various  vortex  trajectories 
will  be  described. 

In  all  cases  described  here,  the  initial  ring  radius  R  was 
used  as  a  representative  length  and  4ttR2/  [  r |  was  used  as  a  time 
scale  to  define  dimensionless  variables.  Let  s  be  a  Lagrangian 
variable  on  the  vortex  loop  which  ranges  from  -1  to  +1  with  s  =  *1 


corresponding  to  the  same  point  (initially  opposite  from  s  =  0  on 
the  ring).  In  dimensionless  form,  equation  (2.20)  becomes 


1  -h 


If  <Vt>  ■  s9n(r)(  jf  +  Otfl  x  - 


(x0-x) 


X0-X|2+p2', 


-  (ff>0  *  (0>oP(s),<iS  +  S9"(r)  [h  R(S-S0)dS 

*  s9"<r>  (|?)0  X  (0)o 


1  ,x  <VX> 

sgn(r)  /_  § f  x  ^ 


-1 


V*l3 


ds 


(3.1) 


image 


For  a  closed  vortex  loop. 


it  may  be  shown 


(Ersoy  and  Walker, 


1985b). 
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where  A  *  X'  •  X'.  The  parameter  u  is  related  to  the  vortex  core 
size  through  equation  (2.17)  and  at  any  instant  the  core  radius 
a(t)  may  be  computed  from  equation  (2.11). 

Starting  from  an  initial  ring  configuration,  the  solution  of 
equation  (3.1)  was  advanced  forward  in  time  using  a  fourth-order 


Runge-Kutta  procedure.  The  spatial  derivatives  were  computed  at  any 
instant  using  central  differences  and  the  numerical  integrations 
were  carried  out  using  Simpson's  rule. 

Some  typical  results  are  shown  in  Figure  3.2  corresponding  to  a 
ring  initially  inclined  at  a  *135°  to  the  plane.  For  this  calcula¬ 
tion,  r  <  0  and  so  the  vortex  loop  recedes  from  the  wall;  the 
initial  core  size  was  aQ  =  .005  and  the  calculation  was  carried  out 
with  a  time  step  At  =  0.005  with  100  points  on  the  vortex.  The  ring 
center  was  initially  at  x  =  0,  y  =  1.2.  As  the  vortex  recedes  from 
the  wall,  there  is  some  contraction  in  the  size  of  the  loop  until 
eventually  the  influence  of  the  wall  becomes  negligible.  Another 
case  using  the  same  numerical  parameters  is  shown  in  Figure  3.3  for 
a  ring  which  recedes  from  the  wall  on  a  normal  trajectory;  note  that 
the  ring  radius  decreases  as  the  influence  of  the  waTl  progressively 
diminishes . 

3.3  Simulation  of  Hydrogen  Bubble  Lines 

The  vortex  loop  trajectories  computed  in  the  previous  section 
are  essentially  unaltered  if  a  uniform  flow  is  superimposed  in  the 
positive  x-direction,  with  the  exception  that  the  vortex  loop  is 
progressively  convected  to  the  right.  In  this  section,  the  main 
interest  is  in  determining  the  distortions  that  will  occur  in  a 
pattern  of  simulated  hydrogen  bubble  lines  near  the  wall  as  the 
vortex  convects  over  the  bubble  lines.  Consider  the  schematic 
diagram  in  Figure  3.4;  here  a  uniform  flow  of  speed  V  is  super- 
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imposed  to  the  right  so  that  the  vortex  loop  is  progressively  con- 
vected  to  the  right.  To  simulate  a  boundary-layer  region  with 
velocity  gradients  near  the  wall  for  y  <  6,  the  unsteady  inviscid 
velocity  distribution  was  multiplied  by  y/<S;  thus  the  velocity  field 
for  y  >  0  is  described  by 


V(x,y,z,t)  = 


—r 

Vinv(x,y,z,t)  y  >  6 


kVinv(x,y,z,t)(|)  y  l  5 


(3.3) 


A  constant  value  of  6  =  0.05  was  used  in  this  simulation.  It  should 
be  noted  that  this  simulation  is  purely  kinematical ;  a  convected 
vortex  will  generally  induce  eruptive  behavior  in  the  viscous  flow 
near  the  wall  (Walker,  1978;  Doligalski  and  Walker,  1984;  Walker 
et  al ,  1987).  The  present  simulation  does  not  take  account  of  this 
effect  which  will  be  addressed  in  Chapter  5. 

Assume  now  that  a  hydrogen  bubble  wire  is  placed  in  the  region 
of  uniform  shear  near  the  wall  at  y  =  hQ  =  6/2.  When  an  electric 
current  is  passed  through  the  wire,  a  line  of  hydrogen  bubbles  is 
produced  in  the  water.  In  the  absence  of  the  vortex,  these  lines 
would  convect  downstream  and  remain  straight  lines;  viewed  from  the 
side  (Figure  3.4(a)).  These  lines  would  appear  as  a  series  of  con- 
vecting  points;  the  top  view  in  the  absence  of  the  vortex  is 
depicted  in  Figure  3.4(b).  In  the  simulations  that  will  be  de¬ 
scribed,  it  was  imagined  that  17  equally  spaced  bubble  lines  were 


placed  in  the  flow  near  the  wall;  the  vortex  was  then  introduced 
into  the  flow  and  the  subsequent  trajectories  of  the  bubble  lines 
were  computed  numerically.  Let  H  denote  the  position  vector  to  a 
point  on  a  simulated  bubble  line;  the  trajectory  of  any  point  on  the 
line  may  be  computed  from. 
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where 


f (y/6 )  = 


y  >  6 


(3.5) 
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y  1  6  . 

In  the  calculations  that  will  be  described  51  mesh  points  were  used 
to  discretize  each  of  the  simulated  bubble  lines.  In  the  figures 
that  will  follow  the  projection  of  the  convecting  vortex  loop  in 
the  xz  plane  will  be  shown  as  a  broken  line  as  indicated  in 
Figure  3.4(b) . 


3.4  Calculated  Results 

Consider  first  the  situation  corresponding  to  a  ring  moving 
away  from  the  wall  on  a  normal  trajectory  (a=n/2)  with  a  uniform 
flow  of  speed  V  =  6  superimposed  to  the  right.  The  trajectory  is 
depicted  in  Figure  3.5  in  side  view.  In  Figure  3.6  the  computed 
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evolution  of  the  simulated  hydrogen  bubble  lines  is  shown  at  various 
stages  of  development.  In  Figures  3.7  and  3.8,  the  development  is 
shown  at  various  stages  as  viewed  from  the  side  and  end  on 
respectively;  in  these  figures  the  dotted  line  represents  the  top  of 
the  simulated  boundary-layer  region  at  y  =  6.  It  may  be  observed, 
upon  inspection  of  Figures  3.5,  3.6  and  3.7,  that  the  motion  induced 
by  the  convected  ring  vortex  acts  to  concentrate  the  simulated 
hydrogen  bubble  lines  into  a  "streaky"  type  of  structure  in  the 
streamwise  direction. 

The  second  case  detailed  here  corresponds  to  the  same  vortex 
configuration  as  in  the  first  case  but  with  the  convection  speed 
reduced  50%  to  V  =  3.  The  vortex  trajectory  in  side  view  is  shown 
in  Figure  3.9.  The  temporal  development  of  the  simulated  hydrogen 
bubble  lines  for  this  case  is  shown  in  Figures  3.10,  3.11  and  3.12 
in  top,  side  and  end  views  respectively.  Again  it  may  be  observed 
that  the  effect  of  the  convected  ring  is  to  draw  the  simulated  flow 
markers  into  a  "streaky"  type  of  structure. 

The  last  case  considered  here  corresponds  to  a  vortex  ring 
which  was  initially  started  inclined  at  an  angle  of  a  =  135c  to  the 
wall.  In  this  case,  the  ring  immediately  starts  to  deform  into  a 
vortex  loop  shape;  the  trajectory  is  depicted  in  Figure  3.13  in  a 
situation  where  the  vortex  is  convected  to  the  right  with  a  speed  of 
V  =  10.  Note  that  in  this  case  the  vortex  stays  closer  to  the  wall 
than  in  the  previous  situations;  as  *  is  increased  toward  180°,  the 
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Side  view  of  the  development  of  simulated 
hydrogen  bubble  lines  due  to  a  convected 

ring  (V= 3) . 
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trajectory  of  the  vortex  loop  remains  progressively  closer  to  the 
wall.  The  deformation  of  the  simulated  hydrogen  bubble  lines  in  top, 
side  and  end  views  is  shown  in  Figures  3.14,  3.15  and  3.16.  Again 
it  may  be  observed  that  the  moving  vortex  collects  the  simulated 
flow  markers  into  a  streaky  structure.  The  results  obtained  in  this 
case  are  representative  of  all  cases  that  were  considered  for 
a  t  90°.  For  trajectories  with  a  closer  to  180°,  the  vortex  remains 
closer  to  the  wall;  the  simulated  bubble  lines  are  stretched  to  an 
increasing  extent  but  also  remain  closer  to  the  wall. 

3.5  Conclusions 

It  may  be  concluded  from  the  calculated  results  of  this  section 
that  a  convected  vortex  ring  or  vortex  loop  is  one  possible  vortex 
motion  which  can  give  rise  to  a  streaky  structure  in  the  flow  near 
the  wall.  In  Chapter  4,  another  type  of  vortex  will  be  considered, 
namely  the  hairpin  vortex,  and  it  will  be  demonstrated  that  this 
type  of  vortex  motion  also  gives  rise  to  a  "streaky"  type  of 
behavior. 


CHAPTER  4 


THE  EVOLUTION  OF  A  HAIRPIN  VORTEX 

4.1  Introduction 

In  recent  times  it  has  become  evident  that  the  hairpin  vortex 
may  be  a  principal  flow  structure  of  the  turbulent  boundary  layer. 
This  viewpoint  was  expressed  first  in  an  early  study  by  Theodorsen 
(1952)  and  a  number  of  contemporary  experimental  studies  have 
supported  this  idea  (see  for  example  Head  and  Bandyopadhyay ,  1981; 
Perry  and  Chong,  1982  and  Acarlar  and  Smith,  1984,  1987a,  1987b). 

In  this  chapter,  the  nature  of  the  evolution  of  a  hairpin  vortex  in 
a  background  flow  is  considered.  The  calculations  described  here 
are  inviscid  simulations  which  are  based  on  an  appropriate  discret¬ 
ization  of  the  Biot-Savart  law.  The  plan  of  the  chapter  is  as 
follows.  In  Section  4.2,  the  equations  of  motion  governing  the 
evolution  of  a  three-dimensional  distortion  in  an  otherwise  two- 
dimensional  vortex  are  described;  the  numerical  scheme  used  to  com¬ 
pute  the  trajectory  of  the  vortex  is  described  in  Section  4.3.  Cal¬ 
culations  for  the  evolution  of  a  haripin  vortex  in  either  a  stagnant 
or  uniform  flow  above  a  plane  wall  are  described  in  Section  4.4. 
Computations  for  the  development  of  a  hairpin  vortex  in  a  shear  flow 
above  a  wall  are  discussed  in  Section  4.5  and  a  discussion  of  re¬ 
sults  is  given  in  Section  4.6.  The  effect  of  moving  vortices  is 
often  observed  experimentally  using  hydrogen  bubble  wires  in  water 
flows;  in  Section  4.7,  the  results  of  a  numerical  simulation  are 
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described  wherein  the  evolution  of  simulated  hydrogen  bubble  lines 
was  computed  numerically  as  a  hairpin  vortex  in  a  shear  flow  con- 
vected  past.  In  Section  4.8,  the  nature  of  the  streamwise  velocity 
fluctuation  induced  by  a  moving  hairpin  vortex  is  described. 

Finally,  calculations  for  the  evolution  of  a  pair  of  interacting 
hairpin  vortices  in  a  shear  flow  are  described  in  Section  4.9. 

4.2  Equation  of  Motion 

It  is  convenient  to  develop  the  equations  of  motion  in  a 
dimensionless  form  and  for  the  hairpin  vortex  motion  of  interest, 
it  is  necessary  to  define  a  velocity  and  length  scale.  Consider 
first  the  situation  depicted  in  Figure  4.1(a)  where  a  two-dimensional 
vortex  of  constant  strength  <  is  located  a  distance  d  above  a  plane 
wall.  For  the  sense  of  circulation  indicated  in  Figure  4.1(a), 
inviscid  theory  predicts  that  the  vortex  will  move  to  the  left  with 
speed, 

■  (4J) 

and  remain  at  constant  height  d  above  the  wall  as  it  convects  in  the 
velocity  field  of  the  image  vortex  below  the  plate.  Now  consider 
the  more  complex  situation  depicted  in  Figures  4.1(b)  and  4.1(c) 
where  a  hairpin  is  in  motion  above  a  plane  wall.  In  the  present 
study,  the  hairpin  was  taken  to  evolve  from  a  small  local  distortion 
in  a  two-dimensional  vortex;  at  large  distances  from  the  hairpin  the 
vortex  will  continue  to  behave  like  a  two-dimensional  vortex,  'or 
the  present  purposes,  it  is  convenient  to  define  dimensionless 
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(a)  Two-dimensional  vortex  (side  view) 
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(b)  Hairpin  vortex  (oblique  view) 


(c)  Hairpin  vortex  (side  view) 

Figure  4.1  -  Vortex  motion  above  a  wall 


variables  in  terms  of  d  and  the  velocity  defined  by  equation  (4.1) 
according  to 


(4.2) 


where  the  asterisks  are  used  to  denote  a  dimensional  quantity.  Note 
that  the  magnitude  of  the  circulation  |rj  is  related  to  the  vortex 
strength  <  by. 


( r  |  =  2tt  k  (4.3) 

For  situations  of  Interest  in  the  present  study,  r  <  0  (c.f.  Figure 
(4.1))  and  the  initial  vortex  configuration  used  is  given  by, 

X(s,t)  =  A{cosai  +  sinaj}e"6s2  +  j  +  sk  .  (4.4) 

This  configuration  corresponds  to  a  two-dimensional  vortex  located 
at  a  distance  1  from  the  wall  with  a  symmetrical  three-dimensional 
distortion  centered  at  s  =  0.  Here  A  represents  the  amplitude  of 
the  distortion  and  a  is  the  angle  the  plane  of  the  distortion  makes 
with  the  plane  wall;  in  addition  g  is  a  (large)  number  whose  value 
determines  the  width  or  spread  of  the  distortion.  The  parameter  s 
is  a  Lagrangian  coordinate  which  ranges  from  to  +°°  along  the 
vortex. 

The  integrals  involved  in  the  Biot-Savart  law  are  over  Infinite 


ranges.  However  it  is  not  possible  to  carry  out  a  numerical  inte- 


gration  for  infinite  values  of  s  and  in  practice,  the  computations 
must  be  carried  out  up  to  some  large  but  finite  value  of  s,  say 
s  =  ±z.  For  l s |  >  it  the  vortex  filament  was  assumed  to  remain  as 
a  straight  line  and 


X(s,t)  =  i  +  j  +  sk  ,  |s|  >  £  .  (4.5) 

Note  that  z  must  be  selected  large  enough  to  ensure  that  there  is 
no  significant  effect  on  the  development  of  the  hairpin  vortex;  as 
the  hairpin  vortex  begins  to  spread  in  the  spanwise  direction,  it 
was  necessary  to  continually  increase  the  value  of  z. 

In  the  present  case,  the  generalization  of  equation  (2.20)  (in 
dimensionless  variables)  for  a  vortex  filament  of  infinite  length  is 

ft  *  X  MV1*  +  Uext  *  {4,6) 

J  * 

where  the  integrals  I.  are  defined  by  the  following  expressions: 

J 

(1)  The  Integrals  on  the  Main  Part  of  the  Vortex 

I-|  =  /  <$v(s,s0,t)ds  ,  T 2  «  jf  ($v(s,s0,t)ds,  (4.7) 
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Here  P  is  defined  by  equation  (2.13),  the  variable  s  is  a  Lagrangian 


coordinate  along  the  vortex  and  the  subscript  o  implies  that  a 


quantity  is  evaluated  at  a  specific  point  s  =  sQ  where  the  local 


velocity  is  being  computed  in  equation  (4.6).  In  addition,  the 


quantity  h  denotes  the  half-interval  removed  from  the  conventional 


integration  over  the  vortex;  for  |s|  <h,  the  representation  of  the 


integrand  is  replaced  by  the  series  expansion  given  in  Appendix  B 


and  there  is 


(2)  The  Integral  Near  the  Field  Point 


I-  =  /  R(s,s  ,t)ds 
J  -h  0 


(4.9) 


where  R  is  defined  by  equation  (2.19)  and  corresponds  to  the  detailed 


series  expansions  near  s  =  sQ  given  in  equation  (B.8)  (Appendix  B); 


(3)  The  Integrals  on  the  Straight  Part  of  the  Vortex 
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where , 
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(4)  The  Integral  of  the  P  function 
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where  P  is  defined  by  equation  (2.13); 


(5)  The  Integral  on  the  Main  Part  of  the  Image  Vortex 


(4.13) 


where. 
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and  finally, 

(6)  The  Integrals  on  the  Straight  Part  of  the  Image  Vortex 


l7  -  /  (Lds  +  /  Q.ds  ,  (4.15) 

-z  z 

where  Q.  is  defined  by  equation  (4.14). 

The  integrals  (4.7),  (4.9)  and  (4.13)  were  evaluated  numerically 
using  a  procedure  based  on  Simpson's  rule;  the  remaining  integrals 
were  evaluated  analytically.  The  integral  appearing  in  equation 
(4.12)  is, 
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where , 


B  =  (M)  .  (1*) 
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(4.17) 


To  evaluate  the  Integrals  in  equation  (4.10),  it  is  first  noted  that 


along  the  straight  portions  of  the  vortex. 


-  (y0-y)i-(x  -x)j 

Q  =  - 5 - 9 -  ,  |s  >  i  ,  (4.18) 

V  {(x-X0)2+(y0-y)2+(z0-s)2+u2}3/2 


where  (xQ,y0  ,zQ)  are  the  coordinates  of  the  field  point  on  the  main 


part  (the  curved  portion)  of  the  vortex  at  s  =  sQ;  x(t)  and  y(t) 


are  the  streamwise  and  normal  coordinates  of  points  on  the  straight 


part  of  the  vortex.  For  fixed  t,  x(t),  y(t)  and  (xQ,y0,z0)  are 
fixed  and  upon  taking  the  integrals  Indicated  in  equation  (4.10)  it 


is  easily  shown  that. 
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In  a  similar  manner,  the  corresponding  integrals  in  equation  (4.15) 


may  be  evaluated  along  the  straight  part  of  the  image  vortex  with 


the  result  that, 
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4.3  Numerical  Methods 


Starting  from  the  initial  vortex  configuration  described  by 
equation  (4.4),  the  subsequent  shape  of  the  vortex  was  tracked  for¬ 
ward  in  time  using  the  standard  fourth-order  Runge-Kutta  algorithm. 
The  initial  vortex  was  split  into  a  number  of  equal  increments  in 
the  variable  s  and  the  subsequent  motion  of  each  point  corresponding 
to  a  fixed  value  of  s  was  tracked  forward  in  time.  The  spatial 
gradients  in  equations  (4.8),  (4.12)  and  (4.14)  were  evaluated  using 
three-point  central  difference  formulae  and  the  integrations  in 
equations  (4.7),  (4.9)  and  (4.10)  were  carried  out  (at  fixed  t) 
using  Simpson's  rule.  After  some  experimentation,  the  value  of  h 
(in  equations  (4.7)  and  (4.9))  used  corresponded  to  one  mesh  length 
in  s. 


4.4  The  Evolution  of  a  Hairpin  Vortex  in  Stagnant  Flow 

In  the  cases  discussed  In  this  section,  the  vortex  was  assumed 
to  have  the  initial  configuration  given  by  equation  (4.4)  and  to  be 
in  an  otherwise  stagnant  flow  above  the  wall  (uext  =  0  in  equation 
(4.6)).  Note  however  that  the  subsequent  development  of  the  vortex 
in  this  case  is  identical  to  that  for  a  vortex  in  a  uniform  flow 
(uext  =  UQi)  exceP*1°n  that  the  vortex  would  be  convected 


to  the  right  (for  UQ  >  0).  In  all  cases  considered  the  vortex  was 


assumed  to  have  a  small  core  radius  a  =  0.02  which  was  held  constant 
during  the  integration.  Note  that  this  may  be  approximation;  the 
theory  of  Moore  and  Saffman  (1972)  applies  only  to  closed  filaments 


and  a  relationship  like  equation  (2.10)  cannot  be  used  since  a 
discontinuity  in  core  size  would  ultimately  develop  at  s  =  ±i.  The 
more  recent  theory  of  Callegari  and  Ting  (1978)  also  applies  only  to 
closed  loop  vortices. 

In  the  first  case  considered,  the  hairpin  vortex  was  assumed  to 
have  an  initial  amplitude  of  A  =  0.5  and  angle  of  inclination  of 
a  =  45°  in  equation  (4.4);  the  value  of  the  spread  parameter  e  used 
was  20.  It  emerges  that  the  self -induced  velocities  predicted  by 
the  Biot-Savart  law  are  rather  large,  particularly  near  s  =  0 
where  the  curvature  of  the  vortex  is  initially  largest.  Consequently 
it  was  necessary  to  take  rather  small  time  steps  in  order  to  avoid 
the  occurrence  of  "wiggles"  or  sharp  corners  in  the  vortex  shape; 
once  sucn  anomalies  appear,  the  numerical  scheme  rapidly  breaks 
down.  A  time  step  of  At  =  0.0002  was  used  in  the  integrations.  For 
the  first  case,  a  uniform  mesh  size  of  as  =  0.005  was  used  and 
initially  the  value  of  sl  employed  was  i  =  2.  Therefore  In  the 
initial  phases  of  the  Integration,  there  were  800  points  along  the 
vortex  from  (-«.,£).  As  t  increased  and  the  hairpin  began  to  spread 
to  the  side,  the  value  of  i  was  increased  to  i  -  3  at  t  *  0.03  and 
then  to  i  =  4  at  t  =  0.06.  The  calculations  are  very  time  consuming 
and  were  stopped  at  t  =  .072  when  the  general  trends  were  evident. 

In  Figures  4.2(a),  4.2(b)  and  4.2(c),  the  evolution  of  the 
hairpin  is  plotted  for  case  1  in  top,  side  and  end  views  respectively; 
the  vortex  is  drawn  at  Intervals  corresponding  to  every  40  time  steps. 


Temporal  development  for  a  hairpin  vortex  in  a 
stagnant  flow.  Case  1  (a=45°,  A=0.5,  6=20); 
the  vortex  position  is  plotted  every  40  time 
steps  (At=0.0002). 
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It  may  be  observed  that  the  head  of  the  vortex  (corresponding  to  the 
tip  of  the  distortion)  moves  rapidly  backward  and  after  a  brief 
interval  starts  to  move  downward  toward  the  wall.  The  rate  of 
movement  of  the  head  should  be  contrasted  with  the  relatively  slow 
backward  movement  of  the  straight  portions  of  the  vortex;  for 
| z j  >  3,  the  vortex  moves  basically  in  the  velocity  field  of  its 
image  and  advances  slowly  upstream.  It  is  evident  that  on  the 
curved  portion  of  the  filament,  the  vortex  evolves  in  such  a  manner 
so  that  high  local  curvature  is  rapidly  diminished.  As  the  head 
bends  back,  the  legs  of  the  vortex  curl  backward  in  a  counter-clock¬ 
wise  direction;  the  disturbance  moves  down  toward  the  wall  and 
spreads  to  both  sides.  With  the  evolution  of  the  legs,  new  second¬ 
ary  hairpin  heads  evolve  to  the  side  of  the  main  distrubance.  The 
calculation  was  terminated  at  t  =  .072  since  the  general  trend  of 
the  disturbance  is  reasonably  well  established  at  this  point;  the 
disturbance  is  expected  to  continue  to  propogate  to  the  side  as  the 
vortex  evolves  into  a  "cork-screw"  shape. 

In  the  second  case  considered,  the  same  geometrical  parameters 
were  used  as  in  Case  1  with  the  exception  that  the  amplitude  of  the 
disturbance  was  doubled  to  A  =  1 .0  and  the  mesh  size  used  was 
as  =  0.01.  The  basic  development  of  the  vortex  is  similar  to  Case  1 
and  is  shown  in  Figures  4.3  at  every  20th  time  step.  It  may  be 
observed  in  Figures  4.3(a)  and  4.3(b)  that  the  head  of  the  vortex 
pulls  rapidly  backward  and  the  vortex  begins  to  evolve  into  the  same 
type  of  "corkscrew"  configuration  observed  in  Case  1.  The  calcula- 
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tion  was  terminated  at  t  =  0.028  since  the  evolution  of  the  filament 
was  judged  to  be  similar  to  Case  1. 

In  the  final  case  considered,  the  effect  of  the  initial  angle 
of  the  disturbance  was  considered  by  calculating  a  case  where  the 
initial  angle  a  =  0°.  The  results  are  shown  in  Figures  4.4.  Again 
the  vortex  head  bends  abruptly  back  and  the  evolution  of  the  vortex 
becomes  similar  to  the  previous  cases.  Consequently  in  all  cases, 
the  initial  distortion  in  the  two-dimensional  vortex  was  observed 
to  evolve  into  a  "corkscrew"  type  of  motion  and  the  disturbance 
spreads  outward  along  the  vortex.  While  the  calculations  were 
carried  out  for  a  vortex  in  a  stagnant  flow,  the  same  vortex 
evolution  would  be  observed  in  a  uniform  flow;  if  a  uniform  flow 
is  superimposed  from  left  to  right,  the  only  difference  in  the  evo¬ 
lution  depicted  in  Figures  4, 2-4. 4  would  be  that  the  vortex  would 
convect  progressively  to  the  right. 


4.5  The  Evolution  of  A  Hairpin  Vortex  in  Shear 

In  the  previous  section,  the  development  of  a  hairpin  vortex 
in  either  a  stagnant  or  uniform  flow  was  considered;  in  this 
section,  the  evolution  of  a  hairpin  vortex  in  a  shear  flow  will 
be  addressed.  The  following  background  flow  for  these  calculations 
was  assumed: 
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Figure  4.4  -  Temporal  development  for  a  hairpin  vortex  in  a 
stagnant  flow.  Case  3  U*0°,  A=0.5,  b=20);  the 
vortex  position  is  plotted  every  30  time  steps 
( At=0.0002).  flQ 


Here  V  is  the  (dimensionless)  uniform  flow  speed  at  large  distances 
from  the  wall  and  y$  is  the  height  of  shear  flow. 

A  series  of  these  calculations  were  carried  out  using  the 
initial  vortex  configuration  given  by  equation  (4.4)  with  A  =  0.5, 
a  =  45°  and  8  =  20;  the  vortex  core  radius  was  a  =  0.02.  The 
evolution  of  the  vortex  was  tracked  by  computing  the  solution  of 
equation  (4.6)  numerically  using  equation  (4.21).  In  all  cases,  the 

dimensionless  height  of  the  shear  flow  was  taken  to  be  y$  =  1.5;  the 

mesh  size  along  the  vortex  was  as  =  0.01.  Initially  the  spanwise 
extent  of  s  was  taken  to  be  at  s  =  £  =  2  and  consequently  there  were 
400  points  on  the  vortex;  as  the  hairpin  vortex  began  to  evolve  in 
the  spanwise  direction  the  value  of  £  was  progressively  increased  to 
£  =  3  and  then  £  =  4. 

The  first  case  for  a  vortex  in  shear  will  be  denoted  here  as 
Case  4  and  for  this  situation,  a  uniform  flow  velocity  (at  large 
distances  from  the  wall)  of  V  =  50  (to  the  right)  was  used;  note  that 

the  velocities  near  the  head  are  about  25  units  to  the  left  (for  the 

initial  configuration  in  equation  (4.4)).  The  value  of  V  selected 
in  this  case  is  thus  about  twice  the  initial  maximum  self-induced 
velocity.  Results  for  Case  4  are  plotted  in  Figures  4.5  every  20 
time  steps;  the  time  step  used  was  At  =  0.0005. 

It  may  be  observed  in  Figures  4.5  that  the  distortion  in  the 
Initial  hairpin  vortex  grows  as  the  vortex  convects  downstream  in  the 
shear  flow.  The  vortex  head  lifts  away  from  the  wall  and  eventually 
rises  into  the  uniform  flow  region.  It  may  also  be  observed  that  the 


vortex  legs  move  progressively  toward  the  wall  and  that  new  hairpin 
disturbances  form  outboard  of  the  primary  disturbance.  It  is 
apparent  that  the  shear  flow  has  the  effect  of  amplifying  the 
initial  distortion  as  well  as  giving  rise  to  additional  hairpin 
vortices  on  the  initially  undisturbed  portions  of  the  vortex.  Cal¬ 
culations  of  this  nature  are  rather  time  consuming  and  were  term¬ 
inated  at  t  =  0.1  for  Case  4  since  the  general  trends  of  the  vortex 
evolution  were  evident  at  this  point. 

In  Case  5,  the  effect  of  Increasing  the  shear  rate  was  studied 
by  taking  V  =  150;  the  background  flow  (given  by  equation  (4.21)) 
is  now  at  higher  speed  and  it  was  necessary  to  decrease  the  time  step 
to  At  =  0.0003  to  obtain  smooth  development  in  the  vortex  trajectory. 
The  results  are  plotted  in  Figures  4.6.  It  may  be  observed  by 
comparison  of  Figures  4.5  and  4.6  that  one  effect  of  the  increased 
level  of  shear  is  to  Increase  the  rate  of  amplification  of  the 
original  disturbance  as  well  as  the  evolution  of  secondary  hairpin 
vortices.  The  trend  for  the  vortex  head  to  lift  out  of  the  shear 
flow  (and  then  bend  back)  Is  evident  in  Figure  4.6(b)  both  for  the 
primary  vortex  head  and  the  secondary  hairpin  vortex  heads.  The 
tendancy  for  the  hairpin  legs  to  move  toward  the  surface  may  also  be 
seen  in  Figures  4.6(b)  and  4.6(c).  It  is  also  clear  from  Figure 
4.6(a)  that  the  increased  level  of  shear  gives  rise  to  an  increased 
level  of  stretching  along  the  vortex  legs. 

The  last  case  considered  corresponds  to  an  even  higher  level  of 
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Temporal  development  for  a  hairpin  vortex  in  a  shear 
flow  (V-150).  Case  5  (a=45°,  A=0.5,  6=20);  the 
vortex  position  Is  plotted  every  20  time  steps 
Ut=0.0003). 
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shear  with  V  =  250  and  is  denoted  as  Case  6;  in  this  case,  it  was 
necessary  to  reduce  the  time  step  to  At  =  0.0002  in  order  to  obtain 
smooth  results.  Calculated  results  for  the  initial  development  are 
depicted  in  Figures  4.7.  It  may  be  observed  that  the  development  is 
similar  to  Case  5  but  that  the  deformation  of  the  vortex  takes  place 
at  a  much  faster  rate.  At  t  =  0.06,  two  secondary  hairpins  have 
developed  outboard  of  the  primary  disturbance,  the  legs  of  the 
vortex  move  progressively  toward  the  wall  (Figure  4.7(c))  and  are 
relatively  close  to  the  wall  at  t  =  0.06.  In  order  to  extend  the 
calculation  to  larger  times,  it  was  necessary  to  extend  the  computa¬ 
tional  domain  from  s  =  Ji  =  3to«.  =  4;a  continuation  of  the  sequence 
for  Case  6  is  shown  in  Figures  4.8.  It  may  be  observed  that  the 
bending  back  of  the  primary  vortex  head  is  now  quite  pronounced.  In 
addition,  the  vortex  legs  have  moved  very  close  to  the  wall;  when  the 
Integration  was  carried  forward  beyond  the  last  stage  depicted  in 
Figures  4.8,  the  vortex  legs  touch  the  wall  after  a  few  time  steps 
and  the  numerical  algorithm  terminated. 

4.6  Discussion 

The  present  results  provide  a  detailed  time  history  of  the 
evolution  of  a  hairpin  vortex  both  in  an  irrotational  flow  as  well 
as  a  shear  flow.  Some  of  the  features  of  the  flow  development  are 
similar  to  those  obtained  by  previous  investigators.  Perry  and 
Chong  (1982),  in  an  attempt  to  discribe  the  mechanisms  of  wall  layer 
turbulence, have  described  the  evolution  of  "A-vortices";  such 
vortices  are  meant  to  be  a  crude  representation  of  hairpin  vortices 
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Figure  4.7  -  Temporal  development  for  a  hairpin  vortex  in  a  shear 
flow  (V=250).  Case  6  (a=45°,  A-0.5,  6=20);  the 
vortex  position  is  plotted  every  30  time  steps 
( At=0 .0002). 
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and  consisted  of  two  finite-length  straight  vortex  filaments 
forming  an  inverted  V  shape  above  the  wall.  A  numerical  integration 
based  on  the  Biot-Savart  law  was  carried  out  for  such  "vortices"  in 
a  stagnant  flow.  Although  this  is  a  rather  crude  and  questionable 
representation  of  a  hairpin  vortex,  the  results  of  Perry  and  Chong 
(1982)  do  show  features  similar  to  the  detailed  integrations 
depicted  in  Figures  4.2  to  4.4;  the  tip  of  the  A -vortex  was 
observed  to  bend  backward  and  the  A  was  observed  to  flatten  and 
spread  in  the  spanwise  direction.  The  present  integrations  are  based 
on  the  full  Biot-Savart  law  with  a  self-consistent  representation  of 
the  vortex  core  flow  and  provide  a  detailed  record  of  the  hairpin 
development. 

Aref  and  Flinchep  (1984)  have  considered  the  evolution  of  a 
hairpin  vortex  in  a  shear  flow  near  a  wall  using  the  localized 
induction  method  due  to  Hama  (1962).  In  this  approach,  the  effect  of 
the  image  vortex  below  the  plate  was  ignored  and  the  self-induced 
velocity  at  a  point  was  computed  from  the  contribution  due  to  a 
small  arc  of  the  vortex  near  the  point;  the  effect  of  portions  of  the 
vortex  remote  from  the  field  point  is  neglected  in  the  "localized 
induction  approximation".  Despite  these  deficiencies,  the  calcula¬ 
ted  results  of  Aref  and  Flinchem  (1984)  do  show  an  evolution  of  a 
hairpin  vortex  in  shear  which  shows  some  qualitative  similarities  to 
the  development  depicted  in  Figures  4.5  to  4.8. 

Recently,  Moin,  Leonard  and  Kim  (1986)  have  computed  the  evo¬ 
lution  of  a  curved  vortex  using  the  Biot-Savart  law;  the  vortex  core 


was  accounted  for  using  an  ad  hoc  model  described  by  Leonard  (1935). 
Their  calculations  were  carried  out  with  relatively  coarse  mesh 
sizes  but  do  show  some  features  that  bear  a  similarity  to  some  of 
the  detailed  results  obtained  here. 

The  present  results  may  be  compared  with  experimental  results 
ofAcarlarand  Smith  (1984,  1981a,  1987b);  in  these  experiments, 
hairpin  vortices  were  created  in  a  sub-critical  laminar  boundary 
layer  and  the  effects  and  motion  of  the  vortices  were  observed  as 
the  hairpins  convected  downstream.  The  vortex  trajectories  that 
have  been  computed  in  this  study  are  in  broad  agreement  with  the 
behavior  observed  in  the  experiments  of  Acarlar  and  Smith  (1984).  The 
heads  of  the  hairpins  were  invariably  observed  to  lift  upward  and 
ultimately  bend  backwards.  The  evolution  of  secondary  hairpin 
vortices  outboard  of  the  original  hairpin  was  clearly  observed  in 
the  experiments  in  agreement  with  the  evolution  shown  in  Figures  4.5- 
4.8.  Figure  4.9  is  reproduced  from  the  study  of  Acarlar  and  Smith 
(1984);  in  this  sequence  the  hairpin  vortex  is  moving  toward  the 
observer  and  as  it  passes  a  hydrogen  bubble  wire  was  pulsed  to 
produce  a  line  of  bubbles  within  the  flow.  In  this  sequence,  the 
movement  of  the  hairpin  legs  toward  the  surface  may  be  observed 
which  Is  in  qualitative  agreement  with  the  evolution  shown  in 
Figures  4. 5-4. 8. 

4.7  Influence  of  a  Hairpin  Vortex  on  Dye  Markers 

A  common  method  of  flow  visualization  is  to  introduce  a 


Figure  4.9  -  Tlow  visualization  sequence  reproduced  from  Acalar 
and  Smith  (19f!4)  showing  a  hairpin  vortex  movinq 
toward  the  observer;  the  visualization  is 
accompl i shed  with  a  hydrogen  bubble  wire. 


sequence  of  hydrogen  bubble  lines  into  the  flow  (in  water)  and  to 
observe  their  subsequent  deformation.  The  wall -layer  streaks  in 
a  turbulent  boundary  layer  are  often  considered  to  be  the  signature 
of  a  moving  hairpin  vortex;  this  view  is  strongly  supported  by  the 
experimental  studies  of  Acarlar  and  Smith  (1984,  1987a,  1987b).  In 
the  present  study,  the  evolution  of  a  hairpin  vortex  in  a  shear  flow 
has  been  calculated  and  in  this  section,  the  influence  of  the  moving 
hairpin  vortex  on  simulated  hydrogen  bubble  lines  in  the  flow  will  be 
eval uated. 

To  carry  out  the  numerical  calculations,  a  total  of  17  equally 
spaced  lines  were  Inserted  in  the  flow  field  at  the  initial  instant 
for  the  vortex  development  of  Case  6  depicted  in  Figure  4.7.  A 
total  of  51  mesh  points  were  used  on  each  half  of  each  line.  In  the 
first  case  considered,  the  lines  were  all  at  an  initial  constant 
height  of  yQ  *  0.8  and  were  separated  by  a  streamwise  distance  of 
0.25  units  in  the  streamwise  direction.  As  the  hairpin  vortex 
evolved  in  the  shear  flow  according  to  the  development  depicted  in 
Figure  4.7,  the  trajectories  of  the  "bubble  lines",  were  computed 
numerically.  A  series  of  trajectories  is  shown  in  Figures  4.10  to 
4.12  in  top,  side  and  end  views;  the  current  vortex  position  is 
shown  as  a  dotted  line.  It  may  be  observed  that  the  simulated 
bubble  lines  are  collected  into  two  nearly  parallel  strips  of 
concentration;  this  behavior  is  closely  consistent  with  the  observed 
results  of  Acarlar  and  Smith  (1984)  (c.f.  Figures  3.21  and  3.25 


run 


of  that  study)  where  horizontal  lines  of  hydrogen  bubble  were 
introduced  at  different  levels  above  the  wall  to  detect  the  typical 
"footprint"  of  a  moving  hairpin  vortex. 

A  similar  sequence  is  shown  in  Figures  4.13  to  4.15  for  the 
"bubble  lines"  located  initially  lower  at  yQ  =0.7.  A  similar 
type  of  development  may  be  observed. 
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(c)  End  view 

Figure  4.15  -  The  simulated  hydrogen  bubble  lines  due  to  a  con- 
vected  hairpin  vortex  (Case  6)  at  t  =  .07;  initial 
bubble  lines  were  located  at  y  =  .7. 


4.8  The  Velocity  Signature  of  a  Hairpin  Vortex 

As  a  moving  hairpin  vortex  convects  in  a  shear  flow,  it 
induces  distortions  in  the  streamwise  velocity  profiles  and  in  this 
section,  the  nature  of  the  induced  flow  will  be  considered  briefly. 
Consider  first  the  situation  corresponding  to  Case  4  of  Section  4.5 
and  select  a  particular  vortex  configuration  at  t  =  0.1;  the  side 
view  of  the  vortex  at  this  stage  of  development  is  shown  in  Figure 
4.16(a).  At  the  eleven  streamwise  locations  indicated  in  Figure 
4.16(a),  the  instantaneous  velocity  fluctuation  profiles  are 
depicted  in  Figure  4.16(b);  these  profiles  represent  the  instantaneous 
velocity  fluctuations  due  to  the  hairpin  vortex  and  were  computed 
using  the  Biot-Savart  law.  It  may  be  observed  in  Figure  4.16(b) 
that  the  velocity  profile  labeled  8  displays  the  largest  negative 
velocity  fluctuation  at  about  y  =  1.5;  as  y  decreases,  the  magnitude 
of  the  fluctuation  decreases.  On  the  other  hand,  the  velocity 
fluctuation  above  y  =  1.5  is  positive  and  peaks  near  1.75.  The 
abrupt  and  relatively  large  reversals  in  velocity  fluctuations  are 
indicative  of  the  presence  of  the  vortex  head.  The  actual  instan¬ 
taneous  velocity  profiles  at  each  station  are  shown  in  Figure  4.16(c). 

Another  perspective  on  the  nature  of  the  streamwise  velocity 
induced  bya  moving  hairpin  vortex  is  shown  in  Figure  4.17  for  Case 
4.  In  this  illustration,  the  instantaneous  velocity  fluctuations 
were  computed  at  a  fixed  streamwise  location  (at  x  =  4.29),  in  the 
time  interval  t  =  0.07  to  t  =  0.12  as  the  vortex  convected  past  the 
measuring  station.  The  vortex  positions  at  various  times  are  shown 


in  Figure  4.17(a)  In  side  view  and  the  instantaneous  velocity 
fluctuations  are  depicted  in  Figure  4.17(b)  at  time  intervals  of 
At  =  0.005.  The  corresponding  velocity  profiles  are  indicated  in 
Figure  4.17(c).  It  may  be  observed  from  Figure  4.17(b)  that  a 
rather  large  velocity  fluctuation  occurs  as  the  head  convects  past 
the  measuring  station.  As  the  hairpin  head  moves  downstream  there 
is  a  negative  velocity  fluctuation  near  the  wall  at  subsequent  times 
which  gradually  diminishes  with  increasing  time. 

A  set  of  instantaneous  velocity  profiles  for  Case  6  of  Section 
4.5  at  fixed  value  of  t  are  shown  in  Figures  4.18;  this  case  cor¬ 
responds  to  the  largest  shear  rate  considered  and  the  instantaneous 
shape  of  the  vortex  in  side  view  is  plotted  in  Figure  4.18(a)  at 
t  =  0.078.  Again  the  abrupt  reversal  in  sign  and  the  relatively 
large  values  of  fluctuating  velocity  are  an  indication  of  the  pre¬ 
sence  of  the  head  at  station  11.  Note  that  in  this  case  the  vortex 
legs  have  moved  relatively  close  to  the  wall  at  this  stage  and  induce 
a  relatively  large  negative  velocity  at  the  wall  (profile  1  in 
Figure  4.18(b)).  In  Figures  4.19,  the  velocity  fluctuation  due  to 
the  moving  hairpin  vortex  are  shown  as  the  vortex  convects  past  a 
fixed  measuring  station  at  x  =  9.2;  the  fluctuations  are  plotted  in 
the  interval  t  =  0.03  to  t  =  0.075  in  Figure  4.19(b).  Again  the 
strong  sign  reversal  in  the  velocity  fluctuation  may  be  observed  as 
the  hairpin  head  passes  the  measuring  station.  As  time  increases, 
a  negative  velocity  fluctuation  is  observed  with  a  peak  that  moves 
gradually  toward  the  wall.  At  t  =  0.065,  the  trailing  portion  of  the 
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Instantaneous  velocity  fluctuations  of  the  measuring  station 
indicated  in  Figure  4.17(a);  note  the  shifted  origins- 
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legs  are  near  the  measuring  station  and  the  negative  fluctuation  is 
at  a  maximum  near  the  wall.  With  the  passage  of  the  hairpin  vortex, 
the  velocity  fluctuation  at  the  wall  decreases. 

The  velocity  distribution  due  to  a  hairpin  vortex  moving  in  a 
shear  flow  close  to  a  wall  has  been  determined  recently  through  use 
of  image  processing  techniques  (Lu,  1987).  In  this  approach,  the 
velocity  profile  and  streamwise  velocity  fluctuations  were  obtained 
at  a  single  measuring  station  by  digitizing  the  high-speed  video 
sequences  of  the  hydrogen  bubble-line  patterns  released  by  a  bubble 
wire  oriented  normal  to  the  wall.  By  comparing  the  results  from 
successive  video  frames,  it  is  possible  to  determine  the  instan¬ 
taneous  streamwise  velocity  distribution  at  the  bubble  wire  (Lu  and 
Smith,  1985).  In  the  experiments  of  Lu  (1987),  this  process  was 
applied  to  hairpin  vortices  shed  from  a  hemispherical  protuberance  in 
a  subcritical  laminar  boundary- layer  flow  on  an  otherwise  smooth 
wall.  Some  results  from  Lu  (1987)  are  shown  in  Figure  4.20  for  the 
streamwise  velocity  and  streamwise  velocity  fluctuations  which  were 
recorded  as  a  hairpin  vortex  passed  by  the  hydrogen  bubble  wire. 

The  sequences  shown  in  Figure  4.20  shown  successive  profiles  with 
time  increasing  from  left  to  right.  The  trends  in  these  results 
are  qualitatively  similar  to  that  shown  in  Figures  4.17(b),  (4.17(c) 
and  4.19(b). 

4.9  Interaction  of  two  Hairpin  Vortices 

In  the  previous  section,  the  motion  of  a  single  hairpin  vortex 


(b)  Velocity  fluctuations 


Figure  4.20  -  Profiles  of  (a)  velocity  and  (b)  velocity 
fluctuation  of  hairpin  vortex  over  a  cycle 
as  obtained  from  hydrogen  bubble  flow 
visualization  pictures  on  plane  of 
symmetry  of  the  hairpin  vortex. 


above  a  plane  wall  was  considered.  In  this  section  the  evolution  of 
two  interacting  hairpin  vortices  is  addressed. 

For  two  hairpin  vortices  above  a  plane  wall,  the  generalization 
of  equation  (4.6)  is. 
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Here  k  =  1 ,2  since  there  are  now  two  hairpin  vortices;  sQ  denotes  a 


point  on  vortex  k  and  the  I.  for  j  =  1 ,7  are  given  by  equations 
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(4.7)  through  (4.20)  for  Integrations  along  vortex  k.  In  addition, 
the  background  flow  ugxt  is  given  by  the  shear  flow  in  equation 
(4.21).  The  two  additional  terms  in  equation  (4.22)  represent  the 
integrations  along  the  second  vortex  and  are  given  by. 


Here  k'  denotes  the  second  vortex  and  the  second  integrand  in 
equations  (4.23)  is  evaluated  along  the  image  of  vortex  k‘. 

The  numerical  method  used  is  similar  to  that  described  for  the 
single  hairpin  vortex.  At  each  stage  in  the  4th  order  Runge-Kutta 
method,  the  intermediate  location  of  each  hairpin  was  computed;  the 
calculation  then  passed  on  to  the  next  sub-step  and  the  process  was 
continued  until  a  complete  time  step  had  been  taken.  Note  that  the 
computing  time  for  two  hairpin  vortices  is  approximately  four  times 
as  great  as  for  a  single  hairpin  vortex.  Because  of  the  relatively 
long  computing  times,  calculations  were  carried  out  for  only  one 
case.  The  two  hairpin  vortices  were  selected  from  two  locations 
computed  in  Case  6  (c.f.  Section  4.5)  with  the  second  vortex 
positioned  slightly  ahead  of  the  first.  The  initial  configuration  is 
shown  in  Figure  4.21;  the  first  vortex  is  labeled  HI  and  is  initially 
located  slightly  behind  the  vortex  labeled  H2.  The  evolution  of 
these  two  vortices  was  then  computed  numerically;  there  were  600 
nodal  points  on  each  vortex  and  the  time  step  used  was  At  =  0.0001. 

The  development  of  the  hairpins  is  shown  in  Figure  4.22  after 
100  time  stepc  at  t  =  0.01.  It  is  evident  that  both  vortices  have 
been  stretched  in  the  shear  flow;  however,  the  legs  of  the  leading 
hairpin  H2  have  moved  closer  together  under  the  action  of  the  hair¬ 
pin  HI  as  well  as  beina  lifted  (Figure  4.22(b)).  The  next  stage  of 
development  is  shown  in  Figure  4.23  where  it  may  be  observed  that 
the  interaction  of  the  two  vortices  has  continued.  The  trailing 


legs  of  H2  have  now  passed  underneath  the  trailing  part  of  HI;  the 
middle  portion  of  the  legs  of  H2  have  passed  upward  between  the  legs 
of  vortex  HI  (Figure  4.23(b)).  The  continuation  of  the  interaction 
may  be  observed  at  the  300th  time  step  at  t  =  0.03  in  Figures  4.24 
and  then  at  t  =  0.4  at  the  400  time  step  in  Figures  4.25;  the  cal¬ 
culations  were  terminated  at  the  480  time  step  and  the  vortex  con¬ 
figurations  are  shown  at  this  stage  in  Figures  4.26.  It  may  be 
observed  from  these  figures  that  the  vortices  tend  to  intertwine  and 
then  reinforce  one  another,  particularly  along  the  vortex  legs.  At 
the  stage  shown  in  Figure  4.26,  the  two  vortices  are  relatively  close 
to  one  another  and  the  calculation  failed  shortly  thereafter.  The 
reason  for  this  is  believed  to  be  due  to  the  fact  that  the  vortices 
were  so  close  that  the  vortex  cores  were  touching.  Once  two  vortex 
cores  approach  this  closely,  it  is  necessary  to  consider  a  local  in¬ 
teraction  region  where  viscous  effects  are  important  and  where  the 
flow  development  is  described  by  the  full  Navier-Stokes  equations. 

Such  a  calculation  has  recently  been  carried  out  by  Weston,  Ting 
and  Liu  (1986)  and  by  Liu,  Ting  and  Weston  (1986)  for  merging  vortex 
rings.  This  type  of  calculation  is  extremely  complicated  and  is 
beyond  the  scope  of  the  present  study.  However  it  is  worthwhile  to 
note  that  when  the  interaction  on  merging  of  vortices  occurs,  the 
vortices  merge  in  interaction  zone  and  then  recombine  and  break  apart. 
The  processes  involved  are  clearly  complex  and  not  well  understood. 

The  interaction  of  a  pair  of  hairpin  vortices  has  been 
observed  experimentally  by  Acarlar  and  Smith  (1984,  1987a,  1987b). 


Hairpin  vortices  which  were  shed  from  a  hemispherical  protuberance 
located  on  the  wall  in  a  subcritical  laminar  boundary  layer  were 
marked  with  dye  and  the  evolution  of  the  vortices  was  observed  as 
the  vortices  convected  downstream.  A  development  similar  to  that 
depicted  in  Figures  4.21-4.26  was  clearly  observed  in  the  experi¬ 
ments.  The  legs  of  the  older  leading  vortex  were  observed  to 
stretch  in  the  shear  flow  but  then  to  move  together  under  the 
influence  of  the  trailing  younger  vortex  and  subsequently  lift 
through  the  legs  of  the  younger  vortex.  The  process  of  vortex 
intertwining  was  clearly  discernible  in  the  experiments  (c.f. 
Figures  3.19  and  3.20  of  Acarlar  and  Smith,  1984). 


CHAPTER  5 


THE  BOUNDARY  LAYER  INDUCED  BY  A  MOVING  HAIRPIN  VORTEX 
5.1  Introduction 

In  Chapter  4,  the  trajectories  of  hairpin  vortices  were  calcu¬ 
lated  for  situations  where  the  hairpin  was  convected  in  both  a  uni¬ 
form  flow  and  a  shear  flow  above  a  wall.  In  this  chapter  the  nature 
of  the  unsteady  viscous  flow  induced  at  the  wall  by  the  moving  hair¬ 
pin  vortex  will  be  considered.  The  motion  of  the  hairpin  vortex 
gives  rise  to  an  inviscid  velocity  distribution  near  the  wall  which 
varies  in  the  streamwise  and  spanwise  direction  and  which  changes 
with  time;  near  the  wall,  a  three-dimensional  unsteady  boundary-layer 
flow  develops  with  time.  The  full  three-dimensional  boundary-layer 
problem  is  complex  and  a  successful  treatment  of  such  problems  would 
require  computational  resources  far  beyond  what  was  available  in  the 
course  of  the  present  investigation.  Consequently  a  method  was 
developed  to  study  the  boundary-layer  development  but  only  on  the 
plane  of  symmetry  of  the  moving  hairpin  vortex. 

The  plan  of  this  chapter  is  as  follows.  In  Section  5.2  the 
nature  of  the  inviscid  velocity  distributions  due  to  the  moving  hair¬ 
pin  vortex  on  the  symmetry  plane  is  examined.  The  boundary-layer 
problem  is  formulated  in  Section  5.3  and  a  set  of  transformations  for 
the  boundary-layer  problem  is  described  in  Section  5.4.  A  discussion 
of  the  numerical  methods  used  is  given  in  Section  5.5  and  calculated 
results  for  one  case  are  described  in  Section  5.6. 


5.2  The  Inviscid  Velocity  Distribution 


In  Chapter  4,  the  three-dimensional  evolution  of  hairpin  vortices 
was  computed  through  a  numerical  integration  of  the  Biot-Savart  law. 

At  any  instant,  the  inviscid  velocity  induced  by  the  hairpin  vortex 
is  given  by, 
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Here  XQ  is  a  vector  to  any  location  within  the  flow  which  is  not 
located  on  the  hairpin  vortex.  The  vector, 


Xh  =  V  +  V  +  zhk  ’  (5-2) 

is  to  an  arbitrary  point  on  the  hairpin  vortex  while  the  vector 


xi  =  xhi  -  yhj  +  zhk  ,  (5.3) 

is  to  the  corresponding  point  on  the  image  vortex.  In  equation  (5.1), 
the  first  and  second  integrals  are  the  contributions  of  the  hairpin 
and  its  image  respectively  to  the  induced  velocity  field  at  XQ. 

Now  consider  the  inviscid  flow  distribution  near  the  wall  which 
corresponds  to  taking  the  limit  y  -*•  0  in  equation  (5.1);  near  the 
wall  the  position  vector  is, 
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Xw  =  xi  +  zk  , 


(5.4) 


and  if  the  inviscid  velocity  distribution  at  the  wall  is  denoted  by. 


uw  =  Uw(x,z,t)i  +  Ww(x,z,t)k  , 


(5.5) 


it  may  readily  be  shown  that, 
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The  form  of  the  velocity  components  near  the  symmetry  plane  is 
obtained  by  taking  the  limit  z  -*■  0  in  equations  (5.6)  and  (5.7);  the 
quantities  x^,  yh  and  az^/as  are  even  functions  of  s  while  zh, 

3xh/as  and  3y^/3s  are  odd  functions  of  s.  Denoting  the  limit  of 
U  as  z  -+  0  by  U  (x,t),  it  is  easily  shown  from  equation  (5.6)  that 
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As  z  ->  0,  Ww  ->  0  but  3Ww/3z  is  nonzero;  defining  eoo(x,t)  by, 


eoo(x,t)  =  lim  -r^-  , 
z-*-0  3 


(5.9) 


so  that. 


it  may  be  shown  using  equation  (5.7)  that, 


ejx.t)  =  -  12sgn(r) 


ax. 

7  *  (>-|1h)ir1  ■  .. 

o  (t*-*^+J'h2«h2,5/Z  h 


(5.11) 


It  is  evident  that  distributions  of  UM(x,t)  and  em(x,t)  may  be 
evaluated  at  any  instant  by  carrying  out  the  indicated  integrations 
in  equations  (5.8)  and  (5.11)  for  several  values  of  x.  This  was 
accomplished  as  follows.  For  a  fixed  value  of  x,  the  integrations  in 
equations  (5.8)  and  (5.11)  were  carried  out  for  |s|  <  i  using  a 
numerical  integration  based  on  Simpson's  rule.  For  |s|  >  a,  the 
vortex  was  assumed  to  be  a  straight  line  filament  and  thus, 
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The  portion  of  the  integral  in  equation  (5.8)  for  s  >  t  may  be  evalu¬ 
ated  by  exact  integration  using. 
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while  the  integrand  for  s  l  in  equation  (5.11)  vanishes. 


(5.13) 


As  examples  of  the  nature  of  the  U  and  6  distributions,  two 
extreme  cases  are  considered  here.  The  first  was  denoted  as  Case  6 
in  Chapter  4  and  corresponds  to  a  hairpin  vortex  embedded  in  a  strong 
shear  flow;  the  temporal  development  of  the  vortex  shape  was  given  in 
Figure  4.7.  In  Figure  5.1,  the  distribution  of  Ujx.t)  is  depicted 
at  various  stages  during  the  evolution  of  the  vortex.  At  t  =  0,  the 
streamwise  velocity  distribution  shows  a  region  of  strong  decelera¬ 
tion  and  then  acceleration  near  the  wall;  in  this  case,  all  values  of 
are  negative.  Note  that  for  the  case  considered  in  Figure  5.1, 
the  basic  shear  flow  is  zero  at  the  wall;  if  the  shear  flow  were  such 
that  the  inviscid  flow  is  reduced  (through  a  region  of  uniform  shear) 
to  a  finite  but  non-zero  velocity  at  the  wall  (say  UQ),  then  the 
only  difference  in  Figure  5.1  would  be  that  all  distributions  are 
relative  to  UQ  rather  than  0. 

As  time  increases,  the  vortex  head  moves  progressively  away 
from  the  wall  and  the  maximum  inviscid  speed  near  the  wall  induced 
by  the  hairpin  diminishes  slightly.  At  a  later  stage  (t  =  0.04), 
the  hairpin  head  has  moved  further  from  the  wall  but  has  also  moved 
further  downstream  as  the  vortex  is  stretched  out  in  the  shear  flow; 
the  influence  of  the  hairpin  head  at  t  =  0.04  may  be  seen  in  Figure 
5.1  where  a  distention  has  developed  in  the  velocity  distribution 
near  x  =  11.  For  increasing  t,  this  distention  amplifies  as  may  be 
observed  at  t  =  0.05  near  x  =  14,  at  t  =  0.06  near  x  =  15  and  then 
near  x  *  18  at  t  -  0.07;  this  effect  is  directly  associated 


Inviscid  streamwise  velocity  on  the  symmetry  plane  near  the  wall 
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with  the  moving  vortex  head  and  this  may  be  confirmed  by  comparing 
Figures  5.1  and  4.7.  At  the  later  stages  depicted  in  Figure  5.1, 
additional  distentions  develop  in  the  inviscid  velocity  distributions 
and  the  origin  of  these  can  be  understood  from  Figure  5.2.  As  the 
vortex  begins  to  develop  secondary  (and  then  tertiary)  heads  which  are 
then  convected  downstream,  local  minima  will  develop  in  the  streamwise 
velocity  distribution;  for  example  the  local  minimum  near  x  =  12  at 
t  =  0.06  (and  near  x  =  13  at  t  =  0.07)  is  associated  with  the  instan¬ 
taneous  location  of  the  secondary  head.  One  final  feature  of  the 
distribution  in  Figure  5.1  is  worthy  of  note;  at  t  =  0.06  and  then  at 
t  =  0.07,  the  maximum  inviscid  speed  has  increased.  This  is  believed 
to  be  due  to  the  influence  of  the  primary  leg  (see  Figure  5.2)  as  it 
moves  closer  to  the  wall  (c.f.  Figures  4.7). 

The  evolution  of  the  spanwise  velocity  is  depicted  in  Figure 
5.3  where  0M(x,t)  is  plotted  at  various  stages  during  the  process. 

It  may  be  observed  that,  as  the  hairpin  is  stretched  in  the  shear  and 
convects  downstream,  there  is  a  spanwise  inflow  toward  the  symmetry 
plane  of  increasing  magnitude  (9oo<0).  Essentially  this  behavior  is 
due  to  the  movement  of  the  legs  toward  the  wall.  It  may  be  observed 
that  the  distribution  is  developing  a  complicated  form  by  t  =  0.06. 
It  is  worthwhile  to  note  that  the  n  distribution  is  not  entirely 

cx< 

negative  and  changes  sign  (even  at  early  times)  at  an  x-location  to 
the  left  of  the  minimum.  One  other  feature  of  interest  is  that  by 
t  =  0.06  the  region  over  which  a  significant  e,  occurs  has  expanded 
in  the  streamwise  direction  considerably  as  the  vortex  stretches  out 
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in  the  shear  flow;  at  this  stage,  there  is  a  significant  in  the 
range  9  <  x  <  18  and  to  some  extent,  the  spanwise  flow  is  similar 
to  that  induced  by  a  pair  of  counter-rotating  streamwise  vortices. 

The  second  case  considered  here  corresponds  to  the  hairpin 
vortex  evolving  in  an  otherwise  stagnant  flow  above  a  plane  wall; 
the  trajectories  of  the  vortex  for  this  case  were  plotted  in  Figure 
4.2.  The  temporal  evolution  of  Uo(x,t)  and  Oc(x,t)  are  plotted  in 
Figures  5.4  and  5.5.  Note  that  the  changes  in  these  distributions 
are  rather  less  dramatic  compared  to  the  situation  when  the  hairpin 
vortex  is  in  shear.  For  a  purely  two-dimensional  vortex,  the 
absolute  minimum  in  U  would  be  at  -4  and  thus  the  effect  of  the 

co 

three-dimensionality  of  the  vortex  is  to  slightly  weaken  the  U^, 
distribution.  If  the  same  hairpin  vortex  were  placed  in  an  otherwise 
uniform  flow  above  the  wall,  the  only  change  from  Figures  5.4  and 
5.5  would  be  that  the  distributions  would  be  convected  to  the  right. 

5.3  The  Boundary-Layer  Problem 

The  boundary  layer  that  develops  near  the  wall  below  the  moving 
hairpin  vortex  is  unsteady  and  fully  three-dimensional.  Let  a  rep¬ 
resentative  length  be  d  (the  distance  of  the  undisturbed  straight 
line  vortex  from  the  wall)  and  take  the  representative  velocity  to  be 
Vc  =  ‘/2d  corresponding  to  self-induced  convection  speed  of  the 
straight  part  of  the  vortex  at  large  distances  from  the  hairpin. 

The  Reynolds  number  for  the  vortex  flow  is  defined  by. 

Re  =  ► /2v  .  (5.14) 


-141- 


Temporal  development  of  the  spanwise  velocity  distribution  near 
the  symmetry  plane  induced  by  a  moving  hairnin  vortex  in  a 
stagnant  flow. 


nrwnnrc  *wv\  .»v\  »vT.  ^vivwwinr 


If  (x,y,z)  are  dimensionless  Cartesian  coordinates  with  corresponding 
velocity  components  (u,v,w),  a  scaled  normal  coordinate  and  velocity 
in  the  boundary  layer  are  defined  by, 

y 1  =  yRe1/2  ,  v‘  =  vRe1/2  ,  (5.15) 
and  the  three-dimensional  unsteady  boundary-layer  equations  become, 
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(5.16) 

(5.17) 

(5.18) 


where  U  and  W  are  the  inviscid  velocities  defined  in  equations  (5.6) 
w  w 

and  (5.7). 

The  numerical  problem  associated  with  the  full  three-dimensional 
boundary-layer  flow  is  very  complex  and  would  require  computer  re¬ 
sources  which  were  not  available  in  the  present  study.  However  the 
boundary-layer  flow  on  the  symmetry  plane  z  =  0  develops  indepen¬ 
dently  of  the  rest  of  the  boundary  layer.  The  flow  near  the  symmetry 
plane  gives  some  insight  into  the  nature  of  flow  near  the  wall  and 
consequently  is  of  interest  here.  Both  u  and  v  are  symmetric  about 
z  =  0  while  w  is  an  odd  function  of  z;  specifically. 


w  %  zo(x,y,t)  + 
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(5.19) 
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where  e  is  a  function  to  be  determined.  Substituting  equation  (5.19) 
into  equations  (5.16)  -  (5.18)  and  taking  the  limit  z  -*•  0,  leads  to 
the  special  form  of  the  boundary-layer  equations  on  the  symmetry 
plane,  viz- 
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(5.22) 


In  normal  circumstances  it  is  convenient  in  a  numerical  solution 
of  unsteady  boundary-layer  equations  to  write  the  streamwise  velocity 
for  example,  in  the  form 


u  =  Uoo(x,t)  F(x,y,t)  .  (5.23) 

The  function  F  is  then  a  normalized  velocity  which  varies  from  zero 
at  the  wall  to  one  at  the  boundary-layer  edge;  the  procedure  is 
successful  when  U  is  such  that  there  are  no  zeros  in  U  for  all  x 

oo  oo 

and  t.  It  was  mentioned  in  Section  5.2,  that  is  negative  for  all 
x  and  t;  this  may  be  confirmed  with  reference  to  Figure  5.1.  The 
distribution  of  0w  however  is  somewhat  more  complex;  it  may  be 
observed  in  Figure  5.3  that  ou  is  predominantly  negative  but  even¬ 
tually  develops  a  region  of  positive  values  (and  two  zeros)  as  the 


legs  of  the  hairpin  approach  the  wall.  It  is  convenient  to  obtain  a 
decomposition  of  of  the  form. 


0oo(x,t)  =  0-|  (x ,t)  +  e2(x,t)  >  (5.24) 

which  is  such  the  and  e2  do  not  have  zeros  for  all  x  and  t. 
Presumably  there  are  a  variety  of  such  decompositions;  however  due  to 
the  complex  nature  of  the  evolution  of  (c.f.  Figure  5.3)  it  proved 
difficult  to  find  a  suitable  decomposition  which  was  valid  for  all  t 
of  interest.  The  following  decomposition  was  obtained  by  trial  and 
error: 
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e2(x,t)  =  ejx.t)  -  G-j ( x ,t) 


(5.26) 


Here  the  subscript  h  denotes  coordinates  along  the  hairpin  and  x^(t) 
is  the  instantaneous  location  of  the  undisturbed  two  dimensional  por¬ 
tion  of  the  vortex  (at  large  distances  from  the  hairpin  disturbance). 
The  temporal  evolution  of  0^  and  02  are  depicted  in  Figure  5.6  when 
it  may  be  observed  that  is  positive  for  all  (x,t)  while  02  is 
always  negative. 

Returning  now  to  the  boundary-layer  equations  on  the  symmetry 
plane,  it  is  convenient  to  introduce  a  functions  <i(x,y,t)  and 
^(x,y,t)  defined  by. 


the  decompositi 


6  =  fy  ’  = 

and 

U  =  3y,V=”fx~<*>  *  ^(x»°»t)=0  .  (5.28) 

The  continuity  equation  (5.22)  is  identically  satisfied  using 
these  variables;  note  however  that  i^(x,y,t)  is  not  a  streamfunction 
and  that  it  is  not  possible  to  define  a  streamfunction  in  a  three- 
dimensional  flow. 


j 


5.4  The  Transformed  Boundary-Layer  Problem 
The  streamwise  coordinate  x  ranges  from 
venient  to  transform  to  a  new  variable  s  defined 


to  °° 
by 


and  it  is  con- 
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Here  xf(t)  represents  the  instantaneous  streamwise  location  of  the 
undisturbed  two-dimensional  portion  of  the  vortex;  this  transforma¬ 
tion  is  one-to-one  and  transforms  the  doubly  infinite  range  in  x  to 
the  range  (-1,1)  in  r,.  The  factors  of  n  in  equation  (5.29)  were 
chosen  for  convenience  with  the  objective  of  transforming  the  stream- 
wise  position  of  the  trailing  legs  of  the  hairpin  vortex  to  near 


C  =  1/2.  Generally  one  objective  of  a  transformation  like  equation 
(5.29)  is  to  expand  or  magnify  regions  where  the  streamwise  velocity 
distribution  is  undergoing  severe  variations.  In  the  present  problem 


this  is  difficult  to  accomplish,  as  may  be  inferred  from  Figure  5.1, 
since  the  locations  of  severe  variation  are  constantly  shifting  in 
time.  The  transformation  in  equation  (5.29)  is  tied  to  the  moving 
two-dimensional  part  of  the  vortex  but  as  the  hairpin  legs  begin  to 
move  near  the  surface,  a  region  of  intense  streamwise  velocity 
develops  at  values  of  x  increasingly  less  than  x  =  x^.  The  trans¬ 
formation  (5.29)  does  a  reasonable  job  of  keeping  this  range  near 
5  =  1/2  in  the  latter  stages  of  the  integrations.  It  was  found  by 
trial  and  error  that  the  form  of  the  transformation  (5.29)  was  very 
important;  furthermore  unless  particular  care  was  taken  to  properly 
expand  the  regions  of  intense  variation  in  Ux  and  6b,  significant 
inaccuracies  occurred  in  the  numerical  boundary-layer  solutions. 

In  the  calculations  that  will  be  reported  here  it  was  assumed 
that  the  effects  of  viscosity  become  important  near  the  wall 
abruptly  at  t  =  0;  for  all  t  >  0,  a  thin  unsteady  boundary-layer  flow 
develops  near  the  wall  as  the  hairpin  vortex  evolves  in  the  shear  flow 
above  the  wall.  In  view  of  this  "impulsive  start"  condition,  it  is 
necessary  to  define  a  new  scaled  normal  coordinate  by 


which  takes  into  account  the  fact  that  the  boundary  layer  initially 
thickens  proportional  to  t  .  In  addition,  the  streamwise  velocity 
and  spanwise  velocity  gradient  in  the  boundary  layer  are  written 
according  to, 


/  J 


u  =  U  —  ,  e  =  e,erfn  +  e,  r1  .  (5.31  ) 
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Here  the  functions  y  and  $  are  related  to  the  original  functions 
and  $  by, 


<p  =  2/t  Uj  ,  $  =  2/t{o.|  I(n)  +  e 2 •:■ }  ,  (5.32) 


and  I(n)  is  the  integral  of  the  error  function  given  by, 


I(n)  =  /  erfndn  =  nerfn  +  —  (e~n  -1)  . 
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(5.33) 


The  vertical  velocity  v  may  be  obtained  by  substitution  of  equations 
(5.32)  into  the  second  of  equations  (5.28)  and, 
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The  functions  v  and  $  are  subject  to  the  boundary  conditions, 


=0  at  n  =  0  , 


for  all  (S,t)  at  the  wall  and  at  the  boundary-layer  edge, 
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for  all  U  .t) . 


Under  these  transformations,  the  streamwise  momentum  equation 
along  the  symmetry  plane  becomes. 
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where , 


U  = 


3n 


and  the  coefficients  in  equation  (5.37)  are  given  by 


and  the  coefficients  in  equation  (5.41  )  are  given  by. 
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Note  that  U^,  6^  and  e2  in  equations  (5.37)  through  (5.43)  as  well 
as  their  derivatives  with  respect  to  5  and  t,  are  to  be  regarded  as 
known  functions  of  z,  and  t  which  may  be  evaluated  from  equations 
(5.8),  (5.25)  and  (5.26);  in  addition,  the  terms  3£,/'<x  and  sc/at  may 
be  evaluated  by  differentiation  of  equation  (5.29),  viz. 
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(5.44  ) 


where  dx;/dt  is  the  speed  (in  the  +  x  direction)  of  the  two-dimen¬ 
sional  part  of  the  vortex.  For  a  given  value  of  z  and  location  x  (t), 
equation  (5.29)  may  be  solved  for  x  and  the  derivatives  in  equation 


(5.44)  may  be  computed. 

The  boundary  conditions  for  equations  (5.37)  and  (5.41)  are 
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\  ab  n 


In  addition,  as  x  +  ±  (£>±1),  the  solution  of  the  boundary-layer 

equations  develops  independently  of  the  solution  for  |x|<“  ( 1 1, i  >  1 ) . 
It  may  be  verified  from  equations  (5.8),  (5.11),  (5.25)  and  (5.26) 
that  U^,  ft-j  and  are  0(x  )  as  |x| it  also  follows  from 

O 

equations  (5.44)  that  the  gradients  h/?x  and  or/ it  are  0(x~  )  as 
|x|-*>j.  Consequently  it  is  easily  confirmed  that  the  limiting  values 
of  the  coefficients  in  equations  (5.39)  and  (5.43)  are. 
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where  A-j(t)  and  B^(t)  are  defined  by, 
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Therefore  if  the  solution  for  U(r,n,t)  and  r (  ■ , n,t)  at  f,  =  *  1 
is  denoted  by  0  and  ",  then  these  functions  satisfy  the 
equations 
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(5.45 


and  the  boundary  conditions  (5.45) 

The  boundary-1 ayer  flow  was  assumed  to  develop  abruptly  from 
rest  corresponding  to  the  insertion  of  the  hairpin  vortex  in  the 
shear  flow  at  t  =  0.  As  t  *  0+,  the  limiting  form  of  equations 
(5.37),  (5.41)  and  (5.43)  is, 

(5.49) 

and  the  solution  of  these  equations  satisfying  conditions  (5.45)  is, 

U  =  T  =  erfn  .  (5.50) 

5.5  Numerical  Methods 

The  numerical  methods  used  in  this  study  are  similar  to  those 
described  by  Ersoy  and  Walker  (1985b).  Starting  from  an  initial 
location,  the  evolution  of  the  hairpin  vortex  was  computed  using  the 
methods  outlined  in  Chapter  4;  the  results  of  this  calculation  were 
stored  on  disk  and  a  boundary-layer  integration  using  a  Crank- 
Nicolson  method  (Ersoy  and  Walker,  1985b)  was  initiated.  The  initial 
boundary-layer  flow  is  described  by  equations  (5.50)  and  to  advance 
the  solution  in  time  the  following  procedure  was  used.  The  solution 
at  ■;  =  *  1  was  first  advanced  one  time  step  by  advancing  the  solution 
of  equations  (5.48);  with  the  solution  at  •,  =  ■  1  known  for  a  given 
value  of  t,  the  solution  of  equations  (5.37)  and  ('  .41)  was  then 


advanced  one  time  step.  Central  difference  approximations  were  used 
for  the  spatial  derivatives  and  the  method  is  second  order  accurate 
both  in  space  and  time.  The  numerical  method  for  the  interior  cal¬ 
culation  (-!<£<!)  is  implicit  and  at  a  given  time  step  iteration  is 
required;  typically  15  to  20  sweeps  of  the  mesh  were  needed  at  each 
time  step  to  obtain  convergence  which  was  deemed  to  have  occurred 
when  two  successive  iterates  for  U  and  t  agreed  to  within  four  sig¬ 
nificant  figures  at  each  mesh  point.  A  typical  time  step  used  was 
At  =  0.0002  and  typical  spatial  mesh  sizes  were  h^  =  0.01  in  the  % 
direction  and  h^  =  0.05  in  the  n  direction.  The  last  of  conditions 
(5.45)  were  applied  at  a  value  of  =  i  =  6  as  an  approximation  and 
this  was  believed  to  be  large  enough  to  ensure  that  there  would  be 
no  significant  change  in  the  solution  with  a  larger  value  of  i.  Con¬ 
sequently  in  a  typical  calculation  there  were  201  mesh  points  in  the 
C  direction  and  121  in  the  direction  normal  to  the  wall;  these  mesh 
sizes  are  believed  to  be  small  enough  to  ensure  good  accuracy. 

The  results  of  the  computations  are  subsequently  depicted 
through  three  types  of  plots,  namely:  (a)  instantaneous  streamline 
plots  in  the  xr,  plane,  (b)  plots  of  lines  of  constant  •  in  the  x- 
plane  and  (c)  limiting  surface  streamlines  in  the  xz  plane  near  the 
symmetry  plane  of  the  hairpin  vortex  at  z  =  0.  The  methods  used  to 
plot  these  results  are  described  in  detail  by  trsoy  and  Walker  (1985b 
but  will  be  discussed  briefly  here. 

As  a  result  of  a  completed  boundary-1 ayer  integration,  values  of 
u  and  v  are  known  at  each  time  step  at  each  point  in  a  rectangular 


mesh  covering  the  xn  plane.  In  general  ,  the  equations  of  an 
instantaneous  streamline  are  given  by. 


dx  _  dy  dz 
u  v  w 


(5.51  ) 


where  T  is  a  parameter  measuring  distance  along  a  given  streamline. 
A  numerical  approximation  to  the  first  of  equations  (5.51)  is 


x 


(5.52) 


where  (uo,vQ)  are  the  'nstantaneous  velocities  at  an  arbitrary 
initial  p..  :nt  (x  ,r,  )  at  fixed  time  t.  By  selecting  a  step  :  :,  the 
streamline  through  the  initial  point  may  be  traced  in  the  xr  plane  at 
time  t;  since  the  spatial  mesh  sizes  h-j  and  h^  are  small  it  is 
possible  to  compute  accurate  values  of  u  and  v  at  any  point  in  the 
plane  using  two-dimensional  linear  interpolation.  In  order  to  ensure 
a  smooth  tracing,  the  step  length  A:  was  selected  so  that  ^Ersoy  and 
Walker,  (1985b), 


(u2+v2)  A-  =  0.005 


(5.53) 


at  any  stage. 

The  lines  of  constant  o  were  obtained  in  a  standard  way.  The 
third  type  of  plot  which  will  be  described  are  limiting  surface 
streamlines  in  the  xz  plane;  at  locations  close  to  the  symmetry 
plane,  the  u  and  w  velocity  components  are  given  by, 


(5.54  ) 


u  =  u(x,y,t)  +  0(z2)  ,  w  =  ze(x,y,t)  +  0(z3), 


where  the  right  sides  of  equations  (5.54)  are  the  symmetry  plane 


quantities  computed  in  this  study.  For  small  z,  equations  (5.54) 


provide  a  reasonable  description  of  the  flow  field  near  the  sym¬ 


metry  plane.  Near  the  surface,  v  is  0(y2)  while  u  and  w  are  0(y) 


thus  for  small  y  the  streamlines  are  described  approximately  by, 


dz  _  dx  _  .T,  „  , 

ze  u  ’  y  ’ 


(5.55) 


where  f  is  a  variable  measuring  distance  along  a  streamline.  In  the 


present  study,  the  limiting  streamlines  were  plotted  using  values  of 


u  and  e  at  y  =  h^  =  O.OSj which  is  one  mesh  length  off  the  surface. 


Note  that  it  is  not  possible  to  clearly  delineate  a  specific  range  of 


values  of  z  for  which  equations  (5.54)  give  a  good  representation 


of  the  velocity  components;  the  maximum  value  of  z  used  here  was 


z  =  0.1.  It  is  also  worthwhile  to  mention  that  an  alternate  repre¬ 


sentation  of  the  limiting  surface  streamlines  is  to  formally  take  the 


limit  as  y  ►  0  in  equation  (5.55)  to  obtain. 


(5.56) 


where 


f(x,t)  = 


>  ' :,y  i  y=0 


u/  ’y I y=o 


(5.57) 


Integration  of  equation  (5.56)  yields. 
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where  xQ  and  zq  are  constants.  In  principle  the  limiting  surface 
streamlines  may  be  plotted  using  equation  (5.58);  in  the  present 
study,  it  was  convenient  to  use  the  same  methods  that  were  used  to 
plot  the  streamlines  on  the  symmetry  plane  to  plot  the  limiting  sur¬ 
face  streamlines  (using  equations  (5.55)). 

5.6  Calculated  Results 

The  computation  of  the  developing  boundary-layer  flow  on  the 
symmetry  plane  is  very  time  consuming  and  each  case  involved  on  the 
order  of  several  days  on  a  VAX  8300.  The  problem  is  complicated  by 
the  fact  that  as  the  hairpin  moves,  the  regions  of  intense  variation 
in  the  and  distributions  constantly  shift  with  time  and  con¬ 
siderable  effort  had  to  be  expended  in  finding  a  suitable  transfor¬ 
mation  of  the  form  (5.29)  which  would  reflect  the  vortex  behavior  in 
a  specific  case.  In  the  calculations  that  will  be  reported>there 
were  201  mesh  points  in  the  £  direction  (h^  =  0.01)  and  121  points 
in  the  n  direction  (h^  =  0.05).  Even  with  the  large  number  of  points 
in  the  £  direction,  it  was  determined  that  unless  the  transformation 
(5.29)  was  adjusted  to  a  specific  case,  that  "wiggles"  would  develop 
in  the  computed  results  as  the  vortex  began  to  severly  deform  in  the 
shear  flow,  these  "wiggles"  were  thought  to  be  a  clear  indication  of 
an  inadequate  streamwise  mesh  spacing  in  particular  regions  of  the 
flow.  Due  to  the  considerable  difficulties  that  were  experienced  in 


calculating  the  boundary-layer  evolution,  a  limited  number  of 
reliable  results  were  obtained  in  the  computational  time  available 
and  results  will  be  reported  here  for  only  one  case.  This  situa¬ 
tion  corresponds  to  Case  6  of  Section  4.5  and  consists  of  a  hairpin 
vortex  convected  in  a  shear  flow  above  a  boundary-layer  flow  near  a 
wall.  The  cases  involving  a  vortex  deforming  in  an  irrotational  flow 
were  not  considered  interesting  and  most  of  the  effort  here  was 
directed  toward  the  situations  involving  vortices  in  shear.  The  time- 
step  used  in  this  calculation  was  At  =  0.0002.  In  the  early  stages 
of  the  integration,  the  difference  equations  were  solved  with 
successive-over-relaxation  iterative  methods;  at  later  times,  it  was 
found  that  under-relaxation  was  more  effective  in  obtaining  a 
converged  solution  of  the  difference  equations  at  each  time  step. 

The  flow  development  at  t  =  0.001  is  depicted  in  Figures  5.7. 

In  these  graphs  (and  all  subsequent  graphs)  the  polygonal  symbol  at 
the  top  of  the  graph  denotes  the  instantaneous  streamwise  location 
of  the  hairpin  head  on  the  symmetry  plane;  the  triangle  denotes  the 
streamwise  location  of  the  undisturbed  two-dimensional  vortex.  In 
Figure  5.7(a)  the  arrows  denote  the  instantaneous  direction  of  flow 
in  the  symmetry  plane  where  it  may  be  observed  that  the  flow  is 
initially  almost  symmetric.  In  Figure  5.7(b)  the  contours  of  con¬ 
stant  e  are  plotted;  positive  values  of  o  mean  that  the  local  span- 
wise  flow  is  away  from  the  symmetry  plane  while  negative  values  of  e 
show  that  the  vortex  is  inducing  a  fairly  strong  inflow  toward  the 
symmetry  plane.  The  limiting  surface  streamlines  in  the  xz  plane 
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Streamlines  in  the  symmetry  plane,  constant  e  contours 
and  surface  streamlines  on  the  wall  at  t  =  0.001. 


are  shown  in  Figure  5.7(c)  where  again,  the  flow  direction  is 
indicated  by  arrows. 

A  further  stage  of  development  is  shown  in  Figures  5.8  at 
t  =  0.01.  It  may  be  observed  that  a  region  of  three-dimensional 
flow  separation  is  now  present;  there  are  two  stagnation  points, 
labelled  SI  and  S2,  associated  with  this  region.  Fluid  particles 
in  the  region  bounded  by  the  limiting  streamlines  through  SI  and  S2 
spiral  toward  a  focus  at  FI  where  the  instantaneous  flow  is  leaving 
the  symmetry  plane.  In  Figure  5.8(b),  it  may  be  observed  that  a 
new  feature  in  the  constant  6  contours  has  developed  in  the  form  of 
a  pocket  of  positive  0  near  the  wall.  The  increasingly  positive 
values  of  e  near  the  center  of  this  region  indicate  increasingly 
strengthening  outflows  from  the  symmetry  plane;  the  center  of  this 
region  is  close  to  the  focus  FI  in  Figure  5.8(a).  The  limiting 
surface  streamlines  at  t  =  0.01  are  shown  in  Figure  5.8(c)  where  it 
may  be  observed  that  limiting  streamlines  through  SI  and  S2  divide 
the  plot  into  three  regions.  In  the  right  of  SI,  the  flow  first 
moves  toward  the  symmetry  plane  but  is  then  deflected  away  toward 
the  limiting  streamline  through  SI.  To  the  left  of  S2  the  flow 
is  toward  the  symmetry  plane  while  between  S2  and  SI  the  surface 
*low  reverses  direction. 

As  time  increases  the  region  of  three-dimensional  separation 
increases  in  extent  in  the  streamwise  direction  as  may  be  seen  in 
Figure  5.9(a)  at  t  =  0.02;  the  stagnation  point  SI  remains  almost 
fixed  relative  to  the  moving  hairpin  vortex  while  the  stagnation 
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Streamlines  in  the  symmetry  plane,  constant  0  contours 
and  surface  streamlines  on  the  wall  at  t  =  0.01. 


point  S2  moves  progressively  downstream  to  the  left.  It  may  also 
be  observed  (in  Figure  5.9(b))  that  the  region  of  positive  outflow 
from  the  symmetry  plane  near  the  wall  is  expanding.  The  expanding 
region  of  recirculating  flow  is  also  apparent  from  the  relative 
position  of  the  stagnation  points  in  the  surface  streamlines  at 
t  =  0.02  in  Figure  5.9(c). 

The  continuation  of  the  trends  observed  in  connection  with 
Figures  5.9  may  be  observed  in  Figures  5.10  at  t  =  0.03.  Note  that 
an  intense  variation  in  the  velocity  field  is  beginning  to  develop 
to  the  right  of  the  region  of  separation  (Figure  5.10(a))  and  also 
that  the  streamlines  to  the  right  of  SI  are  beginning  to  lift  away 
from  the  wall  (in  response  to  the  development  of  the  secondary  hair¬ 
pin  heads  in  the  spanwise  direction).  The  next  stage  of  development 
at  t  =  0.04  is  shown  in  Figures  5.11;  with  the  arrival  of  the  trail¬ 
ing  legs  of  the  hairpin,  the  streamlines  that  were  observed  to  have 
been  lifting  in  Figure  5.10(a)  are  now  seen  to  be  accelerated  in 
Figure  5.11(a)  to  the  right  of  SI.  Note  that  the  region  of  recircu¬ 
lating  flow  has  expanded  to  the  extent  that  the  stagnation  point  S2 
has  now  moved  out  of  the  picture  to  the  left  in  Figures  5.11(a)  and 
5.11(c).  In  addition  the  surface  streamlines  near  SI  are  essentially 
tangential  to  the  symmetry  plane. 

The  flow  development  at  t  =  0.05  is  depicted  in  Figures  5.12; 
the  region  of  recirculation,  depicted  in  Figure  5.12(a),  has  con¬ 
tinued  to  expand  and  has  grown  substantially  in  the  direction 
normal  to  the  wall.  The  intense  variation  in  the  flow  field  to  the 
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Figure  5.9  -  Streamlines  in  the  symmetry  plane,  constan 
surface  streamlines  on  the  wall  at  t  =  0.0 
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Streamlines  in  the  synmetry  plane,  constant 
surface  streamlines  on  the  wall  at  t  =  0.04 


right  of  SI  is  evident  from  the  relative  closeness  of  the  stream¬ 
lines  there.  It  may  also  be  observed  in  Figure  5.12(b)  that  a  pro¬ 
nounced  hump  has  developed  in  the  constant  e  contours  to  the  left  of 
SI;  an  additional  lifting  in  these  contours  may  be  observed  to  the 
right  of  SI . 

The  development  at  t  =  0.06  is  plotted  in  Figures  5.13.  It 
may  be  observed  in  Figure  5.13(a)  that  the  recirculating  region  has 
extended  out  of  the  plot  window  in  a  direction  normal  to  the  wall. 

The  streamlines  to  the  right  of  SI  are  very  close  to  the  wall  and  it 
may  be  observed  from  Figure  5.13(c)  that  at  this  stage  (at  least) 
there  are  no  stagnation  points  to  the  right  of  SI.  It  may  also  be 
seen  that  the  kink  in  the  constant  0  contours  which  was  noted  in 
Figure  5.12(b)  has  now  extended  all  the  way  across  the  boundary 
layer.  Consequently  there  Is  a  narrow  band  of  outflow  from  the 
symmetry  plane  which  extends  all  the  way  to  the  inviscld  region  and 
which  forms  in  response  to  the  development  of  a  small  region  of 
positive  (c.f.  Figure  5.3). 

The  last  stage  of  the  present  calculations  Is  shown  in  Figures 
5.14  at  t  =  0.07.  At  this  point,  there  are  now  two  additional 
stagnation  points  S3  and  S4  which  may  be  observed  in  Figures  5.14(a) 
and  5.14(c).  The  limiting  streamlines  through  S3  and  S4  bound  a 
region  of  strong  outflow  from  the  wall  which  is  fed  by  spanwise  flow 
toward  the  symmetry  plane.  Note  that  an  additional  small  region  of 
positive  0  may  be  observed  to  the  right  of  SI  in  Figure  5.14(b). 

It  may  be  observed  in  Figure  5.14(b)  that  "wiggles"  are  starting 


Streamlines  in  the  symmetry  plane,  constant 
surface  streamlines  on  the  wall  at  t  =  0.05 


Streamlines  in  the  symmetry  plane,  constant  e  contours  and 
surface  streamlines  on  the  wall  at  t  =  0,06. 


to  develop  in  the  constant  e  contours  and  this  is  an  indication  that 
the  numerical  mesh  is  becoming  inadequate.  At  this  stage  the  hair¬ 
pin  vortex  legs  are  rather  close  to  the  wall  (c.f.  Figure  4.8)  and 
the  inviscid  velocity  distributions  are  fairly  complex  (c.f.  Figures 
5.1  and  5.3).  By  t  =  0.07,  the  legs  of  the  vortex  are  moving 
fairly  close  to  the  wall  and  the  inviscid  calculation  is  nearing  the 
end  of  its  range  of  validity.  For  these  reasons,  the  boundary-layer 
calculation  was  not  extended  much  beyond  t  =  0.07.  However  it 
may  be  inferred  from  the  present  results  that  strong  boundary-layer 
growth  is  occuring  as  a  consequence  of  the  imposed  pressure  distri¬ 
butions  due  to  the  moving  hairpin  vortex.  This  point  is  illustrated 
in  Figure  5.15  where  the  temporal  evolution  of  the  displacement 
thickness, 

oo 

6*  =  2/t  /  (l-U)dn  ,  (5.59) 

o 

on  the  symmetry  plane  is  plotted.  The  initial  hump  that  begins  to 
develop  and  grow  in  Figure  5.14  is  due  to  the  expanding  region  of 
recirculating  flow  that  develops  on  the  symmetry  plane  at  a  fairly 
early  stage  in  the  integrations.  The  second  and  third  humps  that 
develop  at  a  later  stage  are  due  to  the  approach  of  the  legs  to  the 
wall  and  the  formation  of  the  secondary  spanwise  hairpins 
respectively. 


In  this  chapter,  a  method  has  been  described  for  the  computa¬ 
tion  of  the  viscous  response  of  a  boundary  layer  to  a  moving  and 
developing  hairpin  vortex  in  a  shear  flow.  It  is  evident  that  the 
moving  hairpin  vortex  gives  rise  to  a  very  complex  unsteady  pressure 
field  which  in  turn  drives  a  complicated  unsteady  boundary-layer 
flow.  The  present  results  describe  only  the  development  of  this 
viscous  flow  on  the  symmetry  plane  but  it  is  evident  that  one  con¬ 
sequence  of  the  evolving  flow  is  that  strong  outward  growth  of  the 
boundary-layer  flow  occurs.  It  is  reasonable  to  expect  that  as  the 
thickening  boundary -layer  flow  penetrates  the  cross  flow  in  the 
inviscid  region  above  the  wall,  a  roll-over  into  a  new  secondary 
hairpin  vortex  will  occur.  This  speculation  is  confirmed  by  the 
experimental  observations  of  Acarlar  and  Smith  (1984,  1987a,  1987b). 


CHAPTER  6 


CONCLUSIONS 

In  the  present  study,  several  aspects  of  the  dynamics  of  hair¬ 
pin  vortices  have  been  examined  from  a  theoretical  standpoint.  The 
results  of  this  study  generally  support  the  suggestions  of  the 
experimental  studies  of  Acarlar  and  Smith  (1984,  1987a,  1987b)  and 
the  idea  that  hairpin  vortices  are  a  basic  building  block  in  turbu¬ 
lent  boundary- layer  flows.  Specific  conclusions  are  as  follows: 

(1)  Hairpin  vortices  are  streak  creators. 

This  is  borne  out  by  the  simulations  described  in  Chater  4 
as  well  as  the  experiments  of  Acarlar  and  Smith. 

(2)  Hairpin  vortices  in  a  shear  flow  develop  secondary  hair¬ 
pin  vortices  which  spread  in  the  spanwise  direction. 

This  was  suggested  by  the  experiments  of  Acarlar  and  Smith  and 
is  confirmed  by  the  calculations  in  Chapter  4;  this  is  a  process  of 
regeneration  wherein  a  hairpin  vortex  can  multiply  itself  in  the 
spanwise  direction. 

(3)  Hairpin  vortices  induce  a  pressure  gradient  which  evokes  an 
eruptive  response  from  the  flow  near  the  wall  and  leads  to  the 
creation  of  secondary  hairpins  through  a  strong  viscous-inviscid 
interaction. 

This  phenomenon  was  clearly  shown  in  the  experiments  of  Acarlar 
and  Smith  in  a  variety  of  situations  and  is  supported  by  the  viscous 
calculations  described  in  Chapter  5.  This  process  is  regenerative 
and  offers  a  plausible  explanation  as  to  how  new  vorticity  from  the 


region  near  the  wall  is  continually  introduced  into  the  outer  layer 
of  a  turbulent  boundary  layer. 

(4)  Pairs  of  hairpin  vortices  act  to  reinforce  one  another 
In  Chapter  4,  it  was  demonstrated  that  two  hairpin  vor¬ 
tices  in  proximity  to  one  another  evolve  in  such  a  way  that  they 
intertwine  and  ultimately  give  the  appearance  of  a  single  hairpin 
vortex.  This  phenomenon  provides  a  possible  explanation  for  the 
phenomenon  of  streak  reinforcement  described  by  Smith  (1984).  lo 
understand  this  mechanism,  assume  that  a  hairpin  vortex  moves  above 
a  wall;  as  it  does  so,  a  streak  occurs  near  the  wall.  After  the  flow 
near  the  wall  is  exposed  to  the  pressure  gradient  induced  by  the 
moving  vortex  for  a  sufficient  period  of  time,  an  eruption  occurs  and 
a  secondary  hairpin  vortex  is  created  In  a  strong  viscous-inviscid 
interaction.  As  the  eruption  occurs  the  streak  is  obliterated  near 
the  point  of  ejection  and  at  this  stage  there  are  two  hairpin  vor¬ 
tices  within  the  flow.  In  view  of  the  present  results,  the  secondary 
and  parent  hairpin  vortex  would  be  expected  to  intertwine  and  rein¬ 
force  one  another.  The  next  observable  phase  of  the  process  should 
be  the  appearance  of  a  new  streak  slightly  downstream  of  the  point 
of  eruption  as  the  pair  of  hairpins  create  a  streak  below  them. 

The  present  study  is  by  no  means  complete  and  there  is  a 
variety  of  research  directions  that  could  be  pursued  in  the  future. 
First,  the  present  vortex  calculation  method  is  believed  to  pro¬ 
duce  very  accurate  results  for  the  motion  of  vortices  having  small 
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cores.  However  the  calculation  method  is  based  on  Moore's  (1972) 
model  of  the  vortex  core  which  is  probably  not  general  enough  to 
describe  all  of  the  vortex  motions  which  are  of  interest  in  connec¬ 
tion  with  the  turbulent  boundary  layer.  The  theoretical  model 
developed  by  Callegari  and  Ting  (1978)  for  closed  vortex  loops  is 
complex  but  is  believed  to  be  on  firmer  ground  than  all  prior  vortex 
calculation  methods;  this  approach  needs  to  be  developed  and  extended 
to  the  types  of  vortices  (of  infinite  extent)  of  interest  in  this 
study.  In  the  second  place,  the  present  calculations  have  shown 
that  hairpin  vortices  evolve  in  complex  ways  which  are  not  easily 
predicted  a  priori  on  an  intuitive  basis;  there  is  therefore  a  need 
to  carry  out  further  numerical  integrations  to  develop  an  improved 
understanding  of  how  hairpin  vortices  react  with  one  another  and 
also  surface  topography  near  walls.  The  present  calculation  method 
is  fairly  time  consuming  however  and  there  is  a  need  to  develop  new 
algorithmic  approaches  which  will  permit  the  accurate  calculation 
the  evolution  of  interacting  vortices  in  three  dimensions  more 
efficiently.  In  the  third  place,  there  is  a  need  to  develop 
approaches  which  will  permit  the  efficient  calculation  and  predic¬ 
tion  of  physical  situations  where  vortex  cores  merge  locally;  in 
the  present  study,  it  was  necessary  to  terminate  the  computation  of 
two  interacting  hairpin  vortices  when  the  two  vortex  cores  touched 
(figure  4.25) . 

A  fourth  area  of  needed  study  is  the  calculation  of  the  full 


unsteady  three-dimensional  boundary  layer  under  a  moving  hairpin 
vortex.  In  the  present  study,  the  boundary- layer  calculations 
were  restricted  to  the  symmetry  plane  of  the  moving  hairpin  vortex; 
these  calculations  do  suggest  that  the  boundary-! ayer  flow  will 
become  eruptive  on  the  symmetry  plane.  However  it  is  entirely 
possible  that  the  boundary-layer  flow  will  evolve  toward  inter¬ 
action  first  at  locations  off  the  symmetry  plane.  In  order  to 
completely  understand  the  viscous  response  at  a  wall  to  a  moving 
hairpin  vortex,  the  full  boundary-layer  problem  must  be  addressed. 
Finally,  the  eruptive  phase  of  the  process,  that  is  observed  in 
experiments  to  lead  to  the  creation  of  secondary  hairpins,  is  of 
considerable  interest.  In  order  to  compute  such  viscous-inviscid 
interactions,  it  will  be  necessary  to  develop  new  algorithmic 
approaches  that  will  permit  the  calculation  of  a  boundary-layer 
flow  interacting  strongly  with  an  outer  inviscid  flow. 
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APPENDIX  A 


In  this  appendix,  it  will  be  demonstrated  that  the  anomolous 
behavior  of  the  integrand  in  Moore's  (1972)  method  which  is  depicted 
in  Figure  2.3  for  the  circular  ring  may  be  associated  with 
inaccuracy  in  the  numerical  differentiation.  For  a  circular  vortex 
ring,  the  position  vector  to  any  point  on  the  ring  may  be  written, 

X(s,t)  =  (Rcosits  ,  Rsinws  ,  Zr)  ,  -1<s<1  ,  (A.l) 

where  R  =  R(t),  Zr  =  Zr(t)  are  the  dimensionless  ring  radius  and 
normal  distance  from  the  wall.  The  first  and  second  derivatives  are 

•> 

—  =  TTf-Rsin-rrs  ,  Rcosws  ,  0}  ,  (A. 2) 

|^r  =  tt 2 { -Rcostrs  ,  -Rsinirs  ,  0}  .  (A. 3) 

If  XQ  denotes  the  position  vector  to  a  point  at  s  =  sQ  on  the  vortex 
where  the  ^df-induced  velocity  of  the  vortex  is  to  be  evaluated, 
then 

XQ  -  X  =  IR(costtso  -  costas)  ,  R(sinTTSQ  -  sinus)  ,  0}  .  (A. 4) 

Suppose  now  that  as  a  result  of  errors  incurred  in  the  numerical 


differentiation,  errors  and  e2  occur  in  the  evaluation  of  aX/as 
and  3 2X l as  which  have  the  following  simple  form, 

-*• 

3  V 

—  =  ( 7t— e-j  )  {-Rsinrrs  ,  RcoSttS,  0}  ,  (A. 5 

32X 

■^7-  =  (7r2-e2)  {-RcOSttS  ,  -Rsimrs  ,  0}  .  (A. 6 

It  is  readily  shown  that, 


(XQ-X)  (  tt- Ei )  R2  ( 1  -cos  it  as) 

3S  (  |X0-X|2+p2}3/2  {2R2(l-cosirAs)+y2}' 


where  as  =  s-s0  and  k  is  a  unit  vector  in  the  z  direction.  Using 
equation  (2.13),  it  may  be  shown  that. 


3X  „  32X  p,  .  .  R2(as)2(.-e1)(»2-e2) 

3S  557  2{R2(4S)2(.-£1)2+U2}3/2 


The  Taylor  series  expansion  of  costtAs  is 


COSttAS  =  1  -  2{  (’"As)2  +  ^-j  (irAs)^  + - , 


and  assuming  that  y2<<ir2R2As2  and  that  e-j  and  e2  are  small,  it  may 
shown  that  the  difference  I  between  equations  (A. 7)  and  (A. 9)  is. 


I  '  TC7T-  ( -3c-,  +  -±)  + 
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It  may  be  observed  from  equation  (A. 10)  that  if  e-| »  ^2  = 


I->0  as  as-*-0  » 


(A. 11) 


which  is  the  required  behavior.  On  the  other  hand,  if  e-j ,  0, 

it  may  be  observed  that  a  large  error  for  small  as  will  occur  in  the 
evaluation  of  the  integrand.  Note  that  for  c^>0  (and  e 2±P)»  the 
sign  of  the  first  term  in  equation  (A. 10)  is  different  from  the 
second  term  and  the  type  of  behavior  depicted  in  Figure  2.3  will 
occur.  It  should  be  remarked  that  this  analysis  is  heuristic  and  a 
complete  error  analysis  was  not  carried  out;  the  analysis  in  this 
Appendix  is  only  an  indication  of  the  potential  problems  that  may 
occur  with  Moore's  (1972)  method. 


-202- 


APPENDIX  B 


In  this  appendix,  the  details  of  the  expansion  of  the  integrand 


R(s,sq)  in  equation  (2.18)  are  described.  If  sQ  denotes  the 


location  on  the  vortex  where  the  self-induced  velocity  is  being 
evaluated  and  s  denotes  any  other  location  on  the  vortex,  a  Taylor 


series  expansion  of  the  position  vector  about  sQ  is. 


X(s.t)  =  x(s0,t)  4  (s-s0)(ff)0  ♦  ^  <$•>„ 


(s-s„)2 


32X\ 


(S-S  )3  3)( 

+  (i_4)  + 

3!  vas3'o 


(B.l) 


where  the  subscript  o  denotes  that  the  derivative  is  evaluated  at 
s  =  sQ;  letting  XQ  denote  X(sQ,t)  it  follows  that 


X  =  -  (s- 


sjr  - 

O  0 


(S-50): 

“7T 


^  - 


U-s0)3 
— J1 - 


X'., 
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where  a  prime  is  used  here  to  denote  differentiation  with  respect  to 
s.  A  similar  expansion  for  aX/ss  is. 


ax 

as 


Xo 


*  <S-S0)X"  * 
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2! 


y  itt 
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(B  .3) 


and  consequently  it  is  easily  shown  that. 


If  x  (X„-X)  =  f  (S-SJ2X;  x  X"  ♦  i  (s-s„)Jx;  X  x; 


/  A  A  A  *  -3  0^1  A  A  A. 
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In  addition. 


!Xq*x|2  =  (s-s 0y%  .  x;  +  (s-s0)3X-  .  X” 
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and  the  binomial  expansion  yields. 
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where , 


a  =  (S-S0)2(X^)2  +  W2  (B .7) 

The  expansion  (B.6)  is  valid  for  (s-sQ)  sufficiently  small. 

An  expression  for  the  Integrand  for  small  (s-sQ)  given  by 
equation  (2.18)  may  now  be  obtained  using  equations  (B.4)  and  (B.6) 
the  result  is. 
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